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^H ABSTRACT 

o 

CN We present a study of dense molecular gas kinematics in seventeen nearby proto- 

^^ stellar systems using single-dish and interferometric molecular line observations. The 

^ non-axisymmetric envelopes around a sample of Class O/I protostars were mapped in the 

N2H+ ( J = 1 ^ 0) tracer with the IRAM 30m, CARMA and PdBI as weh as NH3 (1,1) 
with the VLA. The molecular line emission is used to construct line-center velocity and 
L^' linewidth maps for all sources to examine the kinematic structure in the envelopes on 

^y\ spatial scales from 0.1 pc to ^1000 AU. The direction of the large-scale velocity gradi- 

c? ents from single-dish mapping is within 45° of normal to the outflow axis in more than 

^L{ half the sample. Furthermore, the velocity gradients are often quite substantial, the 

O average being ^2.3 km s~^ pc~^ . The interferometric data often reveal small-scale ve- 

^— ^ locity structure, departing from the more gradual large-scale velocity gradients. In some 

Ci cases, this likely indicates accelerating infall and/or rotational spin- up in the inner enve- 

lope; the median velocity gradient from the interferometric data is ^10.7 km s~^ pc~^ . 
^ In two systems, we detect high-velocity HCO^ ( J = 1 ^ 0) emission inside the highest- 

^-H velocity N2H+ emission. This enables us to study the infall and rotation close to the 

^^ disk and estimate the central object masses. The velocity fields observed on large and 

'^ small-scales are more complex than would be expected from rotation alone, suggesting 

f-^ that complex envelope structure enables other dynamical processes (i.e. infall) to affect 

O the velocity field. 
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Introduction 



Infall and rotation in dense cores and protostellar envelopes both play important roles in the 
formation of protostars and their surrounding disks. Infall must be taking place in the envelopes 
because we observe newborn protostars embedded within their natal clouds. The two classic ana- 
lytic theories of collapse describe infall as either being outside-in ( Larson|[l969 ) or inside-out (Shu 



1977p . The initial angular momentum of the protostellar cloud governs the formation and sizes of 



the proto-planetary disks (e.g. Cassen &; Moosman 1981; Terebey et al. 1984) and will affect the 



ability of the cloud core to fragment into multiple stellar systems (e.g. Burkert &; Bodenheimer 



1993; Bonnell fc Bate||1994 ). The ubiquitous formation of proto-planetary disks ( Haisch et al.|2001 



Hernandez et al. 2007) and the prevalence of binary systems (Raghavan et al. 2010) lends strong 



indirect evidence for the presence of rotation. Therefore, observing these two processes in protostel- 
lar envelopes has enormous potential for constraining star formation theory through comparisons 
to analytic models and numerical simulations. Dense protostellar cores are well known as sites of 



isolated, low-mass star formation (e.g. Shu et al. 1987; Benson & Myers 1989;[McKee fc Qstriker 



2007|) and are the ideal place to study the kinematic structure from large, core/envelope scales 
(~0.1 pc) down to scales near the disk radius. 



Several previous studies have attempted to characterize the rotation in dark clouds and dense 



cores. 



Arquilla &; Goldsmith (1986) examined the kinematic structure of dark clouds using CO on 



40' scales, finding that some exhibit possible rotation signatures derived from velocity gradients. 



Later, Goodman et al. (1993) and Caselli et al. (2002) used the dense molecular tracers NH3 and 



N2H+ to examine rotation in the dense cores within dark clouds (^2-3' scales). Goodman et al. 
(1993) and Caselli et al. (2002) found typical velocity gradients of ^1 to 2 km s~^ 



pc 



which were 



interpreted to be slow, solid-body core rotation. However, the observations had limited resolution 
(60''-90''), nearly the size of the envelopes in some cases. Rotation is not necessarily solid-body; 
however, this assumption simplifies analysis in terms of equipartition and the early data did not 
warrant more sophisticated models. 



Caselli et al. (2002) noted that their finer resolution as compared to Benson & Myers (1989) 



revealed clear deviations from linearity in the velocity gradients. This indicates that the velocity 
structure of cores may be more complex than simply axisymmetric, solid-body rotation. [Chen 



et al. (2007) carried out a higher resolution study using interferometric observations of N2H^ from 
OVRO on a sample of protostars in the Class ( Andre et al.||1993 ) and Class I phases ( |Lada||1987 ). 
They found substantially larger velocity gradients in the inner envelopes, with complex velocity 
fields, only showing probable rotation signatures in a few cases. 

While the possible rotation has been probed using optically thin tracers, the study of infall in 
protostellar envelopes has been limited to observations of optically thick tracers toward the envelope 



center. Zhou et al. (1993) and Myers et al. (1995) used observations of CS and H2CO to infer the 



detection of infall in the envelope of the Class protostar B335. The optically thick molecular lines 
are self absorbed at the line-center and the key signature of infall is the blue-shifted side of the line 



being brighter than the red-shifted side. Later, Di Francesco et al. (2001) found inverse-P Cygni 



hne profiles in H2CO which are interpreted as infall. Furthermore, infall at large-scales may have 
been detected in the pre-stellar core L1544 ( Tafalla et al.|[T998 ). 



These previous studies have interpreted the kinematic structure in terms of axisymmetric 



envelopes. However, our recent studies (i.e. Tobin et al. 2010b; Looney et al. 2007) have revealed 



that the dense envelopes surrounding the youngest, generally Class 0, protostars often have complex, 
non- axisymmetric morphological structure (also see Benson fc Myers||1989 ; Myers et"aL]|1991 ; Stutz 
et al.|[2009 ). Out of the 22 protostellar systems exhibiting extinction at 8/im, only three appeared 
to be roughly axisymmetric. We suggested that the asymmetric envelope structure may play a 
role in the formation of binary systems by making fragmentation easier and that the initial disk 
structure may be perturbed from uneven mass loading. However, in order to understand the effects 
that complex envelope structure has on disk formation and fragmentation, the kinematics of the 
dense gas must be characterized. 

Given the apparent prevalence of complex morphological structure in protostellar envelopes, 
the kinematics of the envelopes must also be studied on scales such that the morphological structure 
is spatially resolved. To examine the kinematic structure of morphologically complex envelopes, 
we have undertaken a molecular line survey focusing on nearby (d < 500 pc), embedded protostars 



drawn from our sample of envelopes in [Tobin et al. ( 2010b| ). We have approached this study 
from two directions. First, we obtained new single-dish N2H+ mapping of sixteen systems with 



the IRAM 30m, a factor of two improvement in resolution over Caselli et al. (2002). Secondly 



we obtained new interferometric N2H^ and NH3 observations of fourteen systems with resolutions 
between 3.5^' and &\ The analysis of both the single-dish and interferometric observations offers 
a comprehensive view of the kinematic and morphological structure from 0.1 pc scales down to 
^1000 AU. Complementary Spitzer imaging gives us a clear view of the outflow angular extent and 
direction in all of these objects. 

In this paper, we are primarily presenting the dataset as a whole and give a basic analysis 
of each object; an upcoming paper will give a more detailed interpretation of the velocity struc- 
tures. The data for each object are quite rich and we find many envelopes with velocity gradients 
normal to the outflow on large-scales in the single-dish maps, persisting down to small-scales in 
the interferometer maps. We have organized the paper as follows: Section 2 discusses the sample, 
observations, data reduction, and analysis. Section 3 presents our general results and discusses each 
object in detail, and Section 4 discusses the basic overall properties of the sample as a whole. 



2. Observations and Data Reduction 

We mapped the protostellar envelopes in the dense gas tracers NH3 and N2H+, which are known 
to be present over a wide range of spatial scales and preferentially trace high density regions (critical 
densities of ^2x10^ cm~^ (Danby et al. 1988; Schoier et al. 2005) and ^1.4x10^ cm~^ (Schoier 



I -4- I 



et al.||2005D respectively), where CO has depleted (e.g. [Bergiii et al.||2002t|Tafalla et al.||2004D . The 



single-dish data were all obtained using the IRAM 30m telescope, mapping N2H+ (J = 1 ^ 0). 
The single-dish observations were followed by interferometric observations of N2H+ ( J = 1 ^ 0) 
using CARMA and the Plateau de Bure Interferometer, in addition to NH3 (1,1) observations using 
the VLA[j including VLA archival data for some objects. We will briefly describe the observations 
and the data analysis procedure for each data set. 



2.1. The Sample 

Our sample of protostellar envelopes selected for kinematic study is directly drawn from the 



objects presented in Tobin et al. (2010b, hereafter Paper I). All protostellar systems observed have 



a surrounding envelope visible in 8/xm extinction. Much of the envelope sample was found within 



archival Spitzer data from the cores2disks (c2d) legacy program (Evans et al. 2009). Given that 
c2d and most other archival observations were short integrations, we were only able to detect very 
dense structures which happened to have a bright backgrounds at 8/xm. Thus, there may be a bias 
towards denser envelopes, but the range of bolometric luminosity is fairly broad with objects <1 
Lq and as much as ^14 Lq. These protostars are listed in Table 1 and are mostly Class systems 
with a few Class Is. The requirement of having an envelope visible in extinction enables us to 
consider the envelope morphology in our interpretation of the kinematic data. We regard the 8/xm 
extinction to more robustly reflect the structure of high density material around the protostars on 
scales > 1000 AU, as molecular tracers may have abundance variations which affect their spatial 
distribution. The requirement of Sfim extinction also biases the sample to more isolated systems, 
only Serpens MMS3 and L673 have multiple neighbors < 0.1 pc away. 



2.2. IRAM 30m Observations 

We observed our sample of protostellar envelopes in two observing runs at the IRAM 30m radio 
telescope on Pico Veleta in the Spanish Sierra Nevada. The first observing run took place between 
2008 December 26 - 30; during this run we only observed four protostars due to poor weather. 
We used the AB receiver system and observed the protostars in the N2H+ ( J = 1 ^ 0) transition 
(z/rest= 93.1737637 GHz; JFiF = 123 -^ 012 component ( |Keto fc Rybicki|[20Tol )); the half power 



beam width (HPBW) of the 30m at this frequency is ^27" . The second observing run took place 
between 2009 October 22-26 using the new Eight Mixer Receiver (EMIR). We observed thirteen 
protostars and revisited L1527 from the first run, using the VESPA auto-correlation spectrometer 
as the backend for all observations. Single-sideband system temperatures of 130 K at v— 93 GHz 
were typical for both observing runs. During the 2008 run we used the 20 MHz bandwidth mode 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 



with 20kHz channels and in 2009 we used 40 MHz with 20kHz channels; see Table 2 the list of 
sources observed and more detail. 

We conducted our observations using frequency-switched on-the-fly (OTF) mapping mode. 
The maps varied in size depending on the extent of the source being observed, most being 3' x 3^ 
Most maps were integrated down to at least ay ^ 150mK for the N2H+ (J = 1 ^ 0) transition, 
noise levels for each map are listed in Table 2. We mapped the sources by scanning in the north- 
south direction and again in the east- west direction to minimize striping in the final map. The scan 
legs were stepped by ^" and we repeated the maps to gain a higher signal-to-noise ratio. Calibration 
scans were taken about every 10 minutes between scan legs and the final maps took approximately 
2 hours to complete. Pointing was checked about every two hours, azimuth and elevation offsets 
were typically ±5'^; the pointing offset remained stable, typically within ^2" during an observation. 
These values agree well with the rms pointing accuracy of ^2" . 

The initial calibration of the OTF data to the antenna temperature scale and CLASS data 
format was performed automatically on-site by the Multichannel Imaging and calibration software 
for Receiver Arrays (MIRA)[j package. Further data reduction was done using CLASS (part of 
GILDAS^). For all molecular lines observed, the frequency switched spectra were folded and baseline 
subtracted using a second order polynomial. We then reconstructed the spectral map on a grid 
such that the FWHM of the beam was spanned by 3 pixels and each pixel is the average of all 
measurements within the FWHM of the beam. 



2.3. CARMA Observations 

The CARMA observations were taken in four observing semesters; most sources were observed 
solely in the D-array configuration (except L1527 and LI 157 see below) which yields ^6^' resolution. 
Three objects were observed in 2009 July and August. The correlator at this time only had three 
bands, operating dual side-band mode, giving 6 spectral windows and was configured for an IF 
frequency of 9L181 GHz. N2H+ (J = 1 ^ 0) and HCO+ (J = 1 ^ 0) (z/rest=89.188518 GHz; 
( Lovas|[r992[ )) were observed in opposite side-bands of one spectral window with 2MHz bandwidth 



and 63 channels giving ^0.1 km/s velocity resolution. The N2H+ ( J = 1 ^ 0) emission spectrum 
comprises 7 hyperfine lines over ^17 km s~^, consisting of two groups of three lines and a single 
isolated line. The 2MHz window is too narrow to observe all the lines; therefore, we observed 
the strongest set of three lines, rather than observing the isolated line to maximize our signal-to- 
noise. The multiple hyperfine transitions enables the optical depth and excitation temperature to 
be determined. The second band was also configured for 2 MHz bandwidth and was centered on 



the main hyperfine component of the HCN ( J = 1 ^ 0) transition (z/rest=88. 6318473 GHz; (Lovas 



1992)), and the third band was configured for continuum observations with 500 MHz bandwidth 



^http://www.iram.fr/IRAMFR/GILDAS 



(1 GHz dual-side band). 

Another set of sources was observed in 2010 April and May. During this time, a new correlator 
was being completed with higher velocity resolution over a wider bandwidth and more spectral 
windows; most data were taken with six bands (twelve spectral windows). For these observations, 
the IF frequency was set to 90.9027 GHz such that N2H+ (J = 1 ^ 0) and HCN (J = 1 ^ 0) 
could be observed in opposite side-bands of the same correlator band. We also observed HCO+ 
(J = 1 ^ 0), H13C0+ (J = 1 ^ 0) (i/rest= 86.754294 GHz), ortho-NHsD (li,i -^ lo,i), (z^rest^ 
85.926263 GHz), and continuum; rest frequencies for these transitions are taken from Lovas| ( [T992 ) . 



All spectral line observations used 8 MHz bandwidth with 385 channels yielding ^ 0.06 km s~^ 
velocity resolution, during reduction this was rebinned to 0.1 km s~^ resolution to reduce noise. 
The continuum observations again had 500 MHz bandwidth (1 GHz dual-side band). See Table 3 
for exact dates of observation for particular sources. 

In addition, L1527 and LI 157 were observed in E-array configuration in 2008 October and 
LI 157 was again observed in D-array in 2009 March. These observations were taken in a three- 
point mosaic pattern to better recover the large-scale emission from the envelopes. The correlator 
was configured with one band for continuum, the other two bands were set to observe N2H+ 
( J = 1 ^ 0) and HCO^ (J = 1 ^ 0). One band had 2MHz bandwidth and was centered on the 
isolated N2H+ component. The other band was configured with 8 MHz bandwidth and 63 channels 
with 0.4 km s~^ resolution to cover all 7 N2H+ (J = 1 ^ 0) hyperfine lines. However, we only 
use the 0.1 km s~^ velocity resolution data for kinematics. The observations of LI 157 are further 



detailed in Chiang et al. (2010). 



All datasets were observed in a standard loop (calibrator-source-calibrator), a bright quasar 
within 15° was used for phase and amplitude calibration. The calibrator was integrated for 3 
minutes while the source was integrated for 15 minutes in each cycle. Absolute flux calibration was 
obtained by observing standard flux calibration sources, see Table 3. Bandpass calibration for the 
continuum bands was accomplished by observing a bright quasar, generally 3C454.3; the spectral 
line bands were bandpass corrected using the noise source. 



Each dataset was processed using the MIRIAD software package ( Sault et al.||1995 ). The raw 



visibilities were corrected for refined antenna baseline solutions and transmission line-length vari- 
ation. The data were then edited to remove uncalibratable data (i.e. poor phase coherence, phase 
jumps, anomalous system temperatures/amplitudes). The bandpass corrections were computed 
using the mfcal routine. The absolute flux calibration was derived using the bootflux routine which 
determines the flux density of the gain calibrator relative to the flux calibration source (absolute 
calibration uncertainty is typically ^10%). The phases and amplitudes were calibrated using the 
mselfcal routine. The phase and amplitude solution calculated for the continuum bands was then 
transferred to the spectral line bands. Continuum images and spectral line cubes were generated 
by inverting the corrected visibilities with natural weighting, creating the dirty map. Then the 
dirty map is CLEANed using the mossdi routine using a clean box of 60^' x 60^' which fits within 



the primary beam of the 10.4m dishes {^72" at A=3.2mm). In this paper, we wih mainly interpret 
the N2H+ data but will comment on the other molecules when relevant and the 3mm continuum 
data are presented in the Appendix. 



2.4. VLA Observations 

The NH3 (1,1) observations were taken with EVLA transition system during the final semester 
of VLA correlator operation in D configuration. Eight-hour tracks were taken for 5 sources dur- 
ing 2009 October, November, and 2010 January (see Table 4). We used the 4 IF mode provid- 
ing two tunings and dual-polarization to observe the NH3 (1,1) and (2,2) inversion transitions 
(z/=23.6944955, 23.7226336 GHz respectively ( |Ho fc Townes||1983D ) with 1.5 MHz bandwidth and 



127 channels yielding ^0.15 km/s velocity resolution. This configuration was able to observe the 
main component and two sets of satellite lines for the NH3 (1,1) transition. We alternately observed 
the source and a calibrator within 15°. Two minutes were spent integrating on the calibrator while 
ten minutes were spent on the target during each cycle. Pointing was updated every hour, 3C84 
was observed as the bandpass calibrator and 3C48 or 3C286 was used for absolute flux calibration. 

The raw visibility data from the VLA were reduced and calibrated using the CASA (Common 
Astronomy Software Applications)|j package. The task importvla was used to convert the VLA 
data to a CASA measurement set. The visibility data were then inspected and edited, specifically 
flagging shadowed data and any VLA antennas. Only EVLA antennas were used in our final dataset 
because we used doppler tracking during the observations, which caused phase jumps between VLA 
and EVLA antennas. Antenna positions were corrected using the gencal task when necessary and 
the absolute flux scale was set using the setjy task (absolute calibration uncertainty is typically 
10%). The phases and amplitudes were calibrated using the gaincal task; the bandpass correction 
was determined using the bandpass task. All the calibrations were then applied using the applycal 
task. We then inspected the corrected data to ensure proper phase correction; any severely outlying 
amplitude points in the source data were also flagged. 

The final spectral datacubes were generated using the clean task. The clean task encompasses 
several individual processes including inverting the visibilities, CLEANing the image, and restoring 
the image. Since the VLA provides a 2' diameter primary beam (field of view), we took several 
additional steps to increase our final image fidelity. We first performed a first-pass CLEANing of 
a spectral data cube over the 2' primary beam. Then we calculated the integrated intensity of the 
main component of the NH3 (1,1) or (2,2) transition by summing the spectral line channels. Next, 
we created a mask image in which pixels below a 2a intensity threshold were rejected and those 
above were kept. This isolated the area of NH3 emission around a particular protostar enabling 
CLEANing down to near a la threshold. Then for the NH3 (1,1) images, we employed the multi- 



^ http://casa.nrao.edu 



scale CLEAN algorithm (Rich et aL]|2008); this algorithm models they sky as the sum of Gaussian 



components of various widths, which for extended sources works better than modeling the sky as 
the sum of many point source CLEAN components. In this paper, we will only present the NH3 
(1,1) data. 

In addition to our own observations, we include VLA archival data[^ for the sources L1527, 
L1521F, and L483. The data for L1521F and L483 were taken in different correlator configurations 
than our observations, see Table 4. The same reduction procedure was applied for these archival 
data as for our observations. 



2.5. Plateau de Bure Interferometer Observations 

LI 157 was observed with the Plateau de Bure Interferometer on 2009 June 17 and 2009 July 8 
in the D-array configuration with 5 antennas operating. It was observed again on 2009 November 
13 in C-array configuration with 6 antennas. During the first track the weather conditions were 
average, with 5 to 10 mm of precipitable water vapor (pwv), respectively and a RMS phase noise 
lower than 64°. During the last two tracks the weather conditions were better, with 3-6 mm and 
5-7 mm of pwv and a RMS noise phase lower than 41° and 24°, respectively 

The 3mm receivers were tuned to the N2H+ (J = 1 ^ 0) transition in the lower-sideband. 
The correlator was configured with two windows for continuum observation, each having 320 MHz 
of bandwidth and one window with 20 MHz bandwidth for the N2H+ (J = 1 ^ 0), yielding a 
velocity resolution of 0.125 km/s. Both the spectral line and continuum observations were taken 
in dual-polarization mode. 

The raw data were calibrated using the CLIC program from the GILDAS software package. A 
standard calibration procedure was used. Visibilities amplitudes and phases for each baseline were 
inspected and bad data (e.g. affected by phase jumps or antenna shadowing) were flagged. Phases 
were then calibrated using observations of the 1927+739 calibrator. Absolute flux was derived from 
observations of MWC349, assuming a flux of 1.15 Jy at 3 mm for that source; uncertainty in the 
absolute flux is ^10%. 

The calibrated data were then reduced using the MAPPING program from GILDAS. A UV 
table was created for the 3mm continuum as well as for the N2H+ (J = 1 ^ 0) line. The 3mm 
continuum visibilities were then subtracted from the line visibilities. Finally, a deconvolved map 
was produced using the the CLEAN algorithm. The synthesized beam in the final map is roughly 
circular with a FWHM of ^ 3.4 x 3.3^' at PA 114°; further details of the observations are listed in 
Table 5. 



^ http://archive.nrao.edu 



2.6. Data Analysis 



2.6.1. Hyper fine Fitting 

The N2H^ and NH3 molecular lines both have a hyperfine emission line spectrum. This allows 
us to robustly determine the line-center velocity and full-width half-maximum (FWHM) linewidths 
in each pixel by fitting all the hyperfine components simultaneously We can fit the line-center 



velocities substantially better than native resolution of the observations. Goodman et al. (1993) 
approximates the observed velocity accuracy as 



CTv 



1.15 



'-peak 



(1.21(5„At;)1/2. 



(1) 



This assumes that the lines have a Gaussian shape, (Jt is the rms noise, Tp^^ak is the peak line 
intensity, Av is the FWHM linewidth, and 6y is the velocity width of the channels. The additional 
factor 1.21 is the intrinsic error in the autocorrelation spectrometer with unity weighting (|Thompson 



et al.||2001 ). Using this relationship for velocity accuracy, if we had a signal-to- noise of 5, assuming 
Av = 0.3 km s~^and 5^=0.1 km s~^, our velocity accuracy would be 0.043 km s~^. Note that this 
is the velocity accuracy for one line, using the multiple hyperfine components of N2H^ and NH3 
we can obtain even greater accuracy. 

To fit the lines, we applied the CLASS hyperfine fitting routines to each line. The hyperfine 
frequencies and line ratios for N2H+ (J = 1 ^ 0) were taken from Keto &; Rybicki (2010) and 



we use the NH3 frequencies and line ratios that are built into CLASS (Rydbeck et al. 1977; Ho &; 



Townes 1983). A semi-automated routine was used for fitting the hyperfine structure across the 



entire spectral map. We first generated an integrated intensity map of the central three hyperfine 
lines, from this map we selected each pixel with 3a detection of N2H^ or NH3 and then generated 
a CLASS script to fit the hyperfine structure at each point and write out a table containing these 
data. The fitting does not require all hyperfine components to have 3a detections, only that 
there are clear detections of the strongest hyperfine lines. Points where the fitting failed due to 
inadequate signal-to-noise were removed from the final table. The table is then used for plotting 
maps of the velocity field, linewidth, optical depth, and excitation temperature, all of which are all 
determined from fitting the hyperfine structure. The zeroth-moment maps (integrated intensity) 
are computed simplistically by measuring the intensity of the hyperfine lines within a given velocity 
range, summing the emission, and multiplying by the channel velocity width. The basic emission 
properties of N2H+and NH3 are separated into three tables: Table[6]for the single-dish N2H+, Table 
[7|for the interferometric N2H+, and Table [s] for the VLA NH3 data. These tables list the line center 
velocities, linewidths, maximum integrated intensities, total optical depth of the transition (sum of 



optical depths for each component), excitation temperature, column density (following [Goldsmith 



& Langer (1999) for N2H+ and (Bourke et al. 1995) for NH3), envelope mass as a function of 



assumed abundance, and approximate envelope radii. 
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2.6.2. Velocity Gradient Fitting 

Using the line center velocities computed from the hyperfine fitting, the velocity gradients are 
computed for each object in three ways. First we simply computed the velocity difference between 
two points offset from the protostar by 10000 AU, normal to the outflow direction, and divided by 
the distance. Second, we computed a linear fits to cuts through the velocity data, taken normal 
to the outflow using points within ±30^' of the protostar and/or within the region of the cloud 
directly associated with the protostar. Lastly, we fit a plane to the entire velocity field using the 



method described in Goodman et al. (1993), but using our own IDL implementation coupled with 



the MPFIT routines ( Markwardt||2009 ). We fit the line-center velocity field with the function 

^LSR = '^0 + <^Aa + bA6 (2) 

where vq is the systemic velocity, Aa and A6 are offsets in right ascension and declination (in arc- 
seconds) and a and b are the velocity gradients per arcseconds in the a and 6 directions respectively. 
The total velocity gradient is then given by 

/ 2 ,2x1/2 206264.98 
^ = (a2 + 62)1/2 >< _ ^3) 

with a position angle (PA) east of north (toward increasing velocity) given by 

9g^tan~^- (4) 

where D is the distance in parsecs and the constant is the number of arcseconds per radian. We use 
all three methods on the single-dish data and only the latter two methods on the interferometer 
data since the extent of N2H+ and NH3 emission varies widely from object to object. 



2.7. Extended Structure Sensitivity 

As a general rule, interferometers filter-out large-scale emission; however, if there is a velocity 
gradient and the lines are well-resolved, the largest scale of emission in a given channel may only be 
a fraction of the full structure. Thus, different portions of the envelopes become visible in different 
velocity channels and the integrated intensity maps (zeroth-moment) build a picture of the envelope 
emission from the multiple velocity components. If all of this emission had resided in one velocity 
channel, substantially more emission would have been filtered-out. Furthermore, the filamentary 
nature of many envelopes enabled them to be viewed over a larger area than an envelope that fills 
the primary beam more evenly. 

We do not combine the single-dish and interferometer data because each probe different size 
scales and it is often advantageous to have the large-scale emission resolved-out of the interferometer 
maps. Resolving-out the large-scale emission isolates the regions of compact emission in the inner 
envelope where there may be significant detail in the velocity field. Furthermore, our analysis of 
only the kinematic structure does not necessitate recovery of all flux; a detailed summary of the 



steps needed to combine single-dish data with CARMA observations is given in Koda et al. (2011) 
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Results 



We mapped the regions around each protostar where we detected the presence of 8/im extinc- 
tion in Paper I. Each envelope was observed to be a bright source of N2H+ ( J = 1 ^ 0) (hereafter 
N2H+) emission in the single-dish data; N2H+ emission was present over much of the area where 
we detect extinction at Sfim in each map, shown in Figures 1 through 23. N2H+ and NH3 (1,1) 
(hereafter NH3) is also detected toward all sources in the interferometer data. The interferome- 
ter observations select out the regions of brightest, compact emission which are usually associated 
with the densest regions of the protostellar envelope. Furthermore, there are many cases where the 
N2H^ or NH3 emission peak is not centered on the protostar in the single-dish and/or interferome- 
ter observations (see Section 4.5). The protostar positions in Figures 1 through 23 are derived from 
their 3mm continuum source, which is often coincident with the 8/xm point source (see Appendix), 
or the 24/xm source where continuum data were not available. 

We know the outflow direction and the angular width of the cavity for all objects in the sample 
from the Spitzer IRAC data. This gives an observational constraint on the region in which the 
outflow may impact the envelope. These data enable the characterization of kinematic properties 
of the envelopes and determination of the origin of the kinematic structure. We can determine 
whether the kinematics reflect the intrinsic velocity structure of the core/envelope or if the outflow 
is likely affecting the observed kinematics. Such distinction is critical to ensure that we are not 
misled in further interpretation. While kinematic information is missing from the Spitzer images 
(i.e. blue and red-shifted sides traced by CO emission), this information is readily available in the 
literature for most objects (Table 1). Note that we can also often infer the blue and red-shifted 
outflow directions from the scattered light morphology and intensity ( Whitney et al.|[2003 ). 



3.1. Similarity of N2H+ and NH3 Emission Properties 

Since our interferometric observations mapped either NH3 or N2H^ for most protostars and our 
single-dish data solely mapped N2H+, it is important to demonstrate that the kinematic structure 
observed in the two tracers is consistent. We have observed the envelope around L1157 in N2H+ 



with the IRAM 30m, CARMA (Chiang et al. 2010), and the PdBI, while observing it in NH3 



with the VLA. The NH3 and N2H+ emission both closely follow the regions of 8/xm extinction as 
shown in Figures [l] and [2] and the velocity maps have a very similar structure. The similar spatial 
emission and kinematic properties indicate that emission from these molecular species arises from 
approximately the same region of the envelope (from 0.1 pc to 1000 AU). Further comparisons 
of our sample can be made to data in the literature: Chen et al. ( |2007 ) for CB230 and IRAS 



03282+3035, L1527 with Goodman et al. (1993), HH211 with Tanner fc Arce (2011), and L483 



with Fuller &; Wootten (2000). In all these cases, the N2H+ and NH3 emission is detected in the 



same regions of the envelope with similar kinematic properties. This conflrms that N2H^ and NH3 
trace similar physical conditions at the level of precision we are probing, in agreement with the 
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results from Johnstone et al. (2010). Section 4.5 further discusses the impact of chemistry on these 
species. 



3.2. Velocity Gradients 

We have computed velocity gradients for all objects using both the single-dish and interferom- 
eter data with the methods described in Section 2.6.2. The velocity gradients calculated from the 
single-dish data are given in Table 9 and the interferometric (N2H^ and NH3) velocity gradients 
are listed in Table 10. The gradient directions from the 2D fitting are plotted in Figures 1-23 and 
also listed in Tables 9 and 10. The one-dimensional (ID) cuts through the velocity fields, normal 
to the outflow and across equatorial plane of each envelope, are shown in Figure [24] presenting an 
alternative view of the envelope velocity structure. Linear fits to the single-dish and interferometer 
data are overlaid on the plots. Notice that in some cases the velocity of the interferometer data 
diverges from the single-dish data. This results from the interferometer filtering-out larger-scale 
emission that dominated the single-dish data and the increased resolution picking out smaller-scale 
velocity structure. 

The gradients calculated for the single-dish data with each method are comparable. The 
median single-dish velocity gradients from the different fitting methods are: 2.1 km s~^ pc~^ (ID 
fitting), 1.7 km s~-^ pc"-*^ (ID two points), and 2.2 km s"-*^ pc"-*^ (2D fitting); the mean gradients 
are 2.3, 2.2, and 2.04 km s~^ pc~^ respectively. The distribution of single-dish velocity gradients 



from the three methods is shown in Figure [25j The lower values of the two point method reflect 
that the region inside the 10000 AU radius of some sources has a higher velocity gradient; there 
are velocity decreases toward the edges of some envelopes that are reflected in the plots in Figure 
24l yielding a preference toward lower gradients in Figure [25| 



The mean velocity gradient of the single-dish sample (^2.2 km s ^ pc ^ ) is about twice the 



average gradient in Goodman et al. (1993) and slightly higher than Caselli et al. (2002), but our 



sample of 16 objects is smaller than their larger samples. We can expect to observe larger velocity 



gradients with our higher resolution data because the lower resolution data in Goodman et al 



( |1993| ) and [Caselli et al. (2002) tend to smear velocity components together. L483 was common 
between our work and the two previous studies, with very similar gradient magnitudes and PAs. 



L1527 and LI 152 were also common between our work and Goodman et al. (1993). The gradient 



directions fit for these sources were similar, but the magnitude of the gradients are different. We 
regard our values as being more reliable because our maps are comprised of substantially more 
independent points. 

The interferometric sample has a median velocity gradient of 8.1 km s~^ pc~^ from 2D fitting 
and 10.7 km s~^ pc~^ from ID fitting, both having a mean gradient of 8.6 km s~^ pc~^ . The 



distribution of interferometric velocity gradients from the two methods is shown in Figure 26 The 
interferometric gradients are often larger than the single-dish gradients by factors of several. The 
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2D fitting method was less reliable for the interferometric data given the often complex velocity 
fields present within the data. Reliable 2D fits could not be obtained for L1157 (PdBI data) or 
Serpens MMS3 due to lack of convergence on their complex velocity fields. Our range of observed 



velocity gradients is comparable to what Chen et al. (2007) found; however, they found three 



protostellar systems with gradients >20 km s ^ pc ^ , we do not find such large gradients in our 
data. 

Using the outflow PAs derived from the Spitzer imaging and CO data from the literature, 
we have compared the gradient PAs (the angle toward increasing velocities) with the outflow PAs 
in Figure [27| The gradient position angles are also marked in Figures 1-23 with solid arrow; the 
majority of the velocity gradients are within 45° of normal to the outflow. The distribution of 
interferometric velocity directions strongly shows a trend for being oriented normal to the outflow 
axis. The gradient directions also generally reflect what is seen at large-scales in the single-dish 



data. Table 10 and Figure 28 show that only three systems have velocity gradient directions which 
differ by more than 45° between the single-dish and interferometric measurements. 

The majority of envelopes in the interferometric sample have an ordered velocity structure. 



despite their often complex morphological structure. In contrast to Volgenau et al. (2006) and 



Chen et al. (2007), many systems have velocity gradients roughly normal to the outflow direction. 



this likely results from our larger sample of observations as compared to the Chen et al. (2007) and 



especially Volgenau et al. ( 2006J ). In addition, visual inspection of the data shows that the large and 
small-scale velocity gradient directions are generally consistent with one another, a feature which 



Volgenau et al. (2006) also sees in their sample. The interferometer observations often reveal small- 
scale kinematic detail near the protostar that is smeared-out in the lower-resolution single-dish 
data. 



3.3. Description of Individual Sources 

We will describe the dataset for each source individually in the following subsections and the 
discuss of our results as a whole is in Section 4. 



3.3A. L1157 



The flattened, filamentary envelope of L1157 has been extensively studied in recent years (e.g. 
Looney et al.||2007 ; Chiang et al.||2010 , Paper I). Its velocity field was first studied in Chiang et al. 
(2010), which showed a weak velocity gradient along the filament, normal to the outflow. We 



subsequently observed L1157 with the IRAM 30m, the PdBI, and the VLA. Figure [T] shows the 
data from the IRAM 30m and PdBI and Figure [2] shows the data from the VLA and CARMA. The 
30m data detects a large-scale velocity gradient normal to the outflow and this gradient follows the 
long axis of the envelope. On the east side of the envelope, where the emission curves downward. 
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the N2H^ continues to trace dense material as it becomes more blue- shifted. 

The PdBI data are shown in the bottom panels of Figure [ij The N2H+ emission appears 
double-peaked on small-scales. The 3mm dust continuum shows that the protostar resides between 
the peaks (see Appendix); the peaks are at radii of ^1000 AU (^3.5^^). The velocity field on scales 
>15'^ reflects what has been seen with the CARMA, 30m, and VLA data. On the other hand, the 
velocity of the gas becomes highly red-shifted (^1 km s~^) at the N2H+ peaks as compared to the 
surrounding gas at ^1000 AU from the protostar. The high- velocity gas was observed in [Chiang 



et al. (2010); however, its spatial location was not well-resolved in the CARMA data due to having 
a factor of ^2 lower resolution than the PdBI data. Note that the 2D velocity gradient fits for the 
VLA and CARMA data have position angles that differ by 80°. This results from the north-south 
gradient being more prominent in the CARMA N2H+ data while the east-west gradient appears 
more prominent in the NH3 data. A 2D gradient could not be fit to the PdBI data due to the 
highly complex velocity field present on small-scales. 

The single-dish velocity and linewidth maps show a gradual large-scale velocity gradient in 
LI 157 with broad linewidth in the inner envelope, consistent with the PdBI linewidth map in 
Figure [ij The broad linewidths were also shown in lower resolution data from Chiang et al. (2010). 



Both the east and west peaks in the PdBI maps have velocity wings toward the red and are not 



significantly extended toward blue-shifted velocities. Chiang et al. (2010) attributed the broad 
inner envelope line wings to infall. However, close examination of our higher-resolution PdBI data 
indicate that the blue and red-shifted line wings may result from outflow interaction effects. Figure 
[3] shows that the most red-shifted emission is shghtly shifted to the southeast, along the outflow 
and traces one edge of the northern outflow cavity. Furthermore, the blue-shifted emission also 
seems to outline the southern outflow cavity quite well. Thus, it appears that the outflow is may be 
entraining inner envelope material, while at larger scales the velocity structure appears unaffected 
by the outflow, tracing the intrinsic kinematic structure of the envelope. The spatial overlap of 
red and blue-shifted N2H+ southeast of the protostar on the blue-shifted side of the outflow can 
understood if the outflow is entraining material within a symmetric cavity in the inner envelope, 
producing both a blue and red-shifted component. There also appear to be outflow effects in the 
N2H+ emission on the northwest side of the protostar but to a lesser extent. This is the first 
example of possible outflow entrainment in such a dense gas tracer. However, we note that rather 
than entrainment, the broad N2H^ emission could also result from an outflow shock at that location. 



3.3.2. LI 165 

The protostar L1165IRS is located within a ^1.5 pc (17^) long filamentary dark cloud (Paper 
I). The narrow filament, from which the protostar has formed, is normal to the protostellar outflow. 
We mapped a 4' section of the dark cloud surrounding the protostar, as shown in the left panels 
of Figure [4J The N2H+ emission is highly peaked very near the protostar, it is slightly offset 
4.87^^(^1462 AU) to the southeast, and there is low-level extended emission associated with regions 
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of 8/iiii extinction. The velocity field from the single-dish N2H^ data shows a fairly linear gradient 
nearly normal to the outflow; in areas away from the protostar the filament generally seems to have 
a fairly constant velocity with little variation and small linewidth. 

The N2H+ data from CARMA reveal a small-scale structure that is extended in the direction 
of the filament axis and the emission is strongly correlated with the small-scale S/xm extinction 
shown in Figure |4j At the edges of detected emission, the velocity field of the interferometer map 
is consistent with the single-dish data. However, near the protostar there is a 0.35 km/s velocity 
shift between the blue and red-shifted velocity peaks that are on opposite sides of the protostar. 
We note that this velocity gradient is not perfectly normal to the outflow, but rather offset by 
about ^30°. However, the most red-shifted emission is slightly extended in the direction of the 
outflow, but the linewidth peak is extended normal to the outflow and the red-shifted emission at 
this location does not appear to be outflow affected. Thus, the gradient from the envelope itself 
appears to be normal to the outflow and tracing gravitationally dominated motion. 

We show the HCO^ emission in Figure [5| in which the red and blue-shifted emission is confined 
to two clumps, that are oriented normal to the outflow and offset from the protostar by 3'^ The 
position- velocity plot of these data show the blue and red-shifted emission extending ~2 and ^1.5 
km s~^ away from the systemic velocity respectively. Given the orientation of the blue and red- 
shifted emission normal to the outflow and a lack of emission along the outflow, the HCO^ emission 
appears to be originating from the inner envelope. Assuming that the HCO^ emission indicates 
rotationally supported motion, we can calculate the enclosed mass using M — R^v^ jG which gives 
Menc ^ 2.OM0. This is not unreasonable for this source which has L^^/ ^14L0 assuming a distance 
of 300 pc. However, if the velocities result from equal contributions of rotation and infall, then 
the enclosed mass would be a more modest 0.5 Mq. We have overlaid lines representing Keplerian 
rotation (or infall) on Figure Isj the O.5M0 curve matches the data much better than the 2.OM0 
curve. 



3.3.3. CB230 
CB230 is an isolated protostar that formed at one end of its natal globule (Paper I). This 



protostar was discovered to be a wide binary system by Yun (1996) with a separation of ^10'^ 
the companion is evident in the 8/xm images shown in Figure [6j The envelope around CB230 was 
classified as a "one-sided" envelope in Paper I due to the 60'^ (30000 AU) extension of the extinction 
envelope to the west, while the 8/xm extinction terminates just ^25'^ (6000 AU) on the eastern side. 
The single-dish N2H^ emission is consistent with the 8/xm extinction observations, the emission 
falling off steeply to the east and more extended to the west. The large-scale extension of material 
beyond the region of detected N2H+ emission suggests that it has lower density and has not formed 
N2H+ at detectable levels, see Section 4.5 for further discussion. 

The NH3 emission from the VLA is strongest in a "bar," about 10° from normal to the outflow. 
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Notably, at the location of the protostar there is a "hole" in the ammonia emission; a similar 
depression of emission is seen in N2H+ ( J = 1 



0) by Chen et al. (2007) and Launhardt et al. 
( |2001| ). The region of decreased emission is ^2200AU in diameter, similar in size to the double- 
peaked N2H+ emission in L1157. The lower level emission on the eastern side of the protostar 
extends northward along the outflow cavity. Incidentally, the northward extension is also where 



the scattered light emission in the near-IR and Spitzer 3.6/xm is brightest (Paper I; [Launhardt et aL 



2010). 



The velocity field of the single-dish N2H+ shows a fairly linear gradient, normal to the outflow, 
with a "plateau" of the most highly red/blue-shifted emission ±30^' from the protostar. South of 
the protostar, the N2H+ linewidth peaks, similar to LI 157. The NH3 velocity gradient from the 
VLA data is in the same direction as the single-dish gradient and similar to the interferometric 



N2H+ map from Chen et al. (2007). While the single-dish gradient is fairly gradual, the gradient 



from the VLA observations has an abrupt shift from blue to red-shifted emission coincident with 
the protostar. The VLA NH3 data also show a velocity "plateau" in the blue and red-shifted 
emission, with the highest relative-velocity emission being ±15'^ from the protostar. The linewidth 
remains fairly constant throughout the regions near the protostar, peaking at 0.5 km/s. The region 
of highest linewidth also corresponds to the region of strongest ammonia emission. The large 
linewidth seen near the protostar in the single-dish data is not reflected in these NH3 data nor the 



N2H+ data of Chen et al. (2007). The line-center velocity changes quite rapidly at the location of 



the large single-dish linewidth; therefore, the linewidth peak is likely due to the unresolved velocity 
gradient. This means that the outflow is not likely to be affecting the kinematics of the inner 
envelope in this protostar. 



3.34. HH108IRS 

The protostar HH108IRS, the driving source of HH108, is located within a large-scale filament, 
■^0.5 pc in length, 1.75° south of the Serpens star forming region (Harvey et al. 20061). There are 



at least two protostars forming in the filament: the higher luminosity object HH108IRS and the 
deeply embedded source HH108MMS ( |Chini et al]|200l[ ). The single-dish N2H+ map in Figure [t] 
shows an emission peak coincident with HH108IRS,but slightly offset from the protostar ^5'^ (1500 
AU). The N2H+ map from CARMA reveals that the N2H+ peak emission is truly offset from the 
protostar and emission is extended normal to the outflow direction, forming a flattened structure 
across the protostar. 

The single-dish N2H^ velocity map indicates that there may be a slight gradient normal to the 
outflow. The CARMA N2H+ velocity map reveals that there is indeed a velocity gradient normal 
to the outflow, though its structure is complex. Southeast of the protostar the velocities are red- 
shifted and moving toward the protostar the velocities are becoming more blue-shifted. This trend 
continues after crossing the protostar and moving northwest, but then the trend reverses itself 
rapidly and becomes more red-shifted. We also note that the linewidth is ^0.6 km/s within ±10'^ 
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(3000 AU) around the protostar, indicative of a dynamic environment near the protostar. The 
single-dish hnewidth is consistent with this value as well. 



3.3.5. HH108MMS 

HH108MMS is the nearby neighbor to HH108IRS, separated by -60'' (-0.09 pc). This proto- 
star is deeply embedded and invisible at 24/im, only becoming visible at 70/xm (Paper I). There are 
no IRAC data available for this object, therefore we are showing an ISPI Ks-band (2.15/xm) image 
from Paper I. Despite the lack of 8/im extinction observations, we can clearly see the dense material 
of the envelope blocking out the rich background star field. The single-dish N2H+ in Figure [s] shows 
a slight extension toward the location of the protostar (derived from 70/xm and 3mm continuum), 
while the CARMA N2H+ observations clearly show the N2H+ emission centrally peaked on the 
protostar. 

The single-dish velocity field does not show much structure, the region within ±30'' (0.045 pc) 
of the protostar has a roughly constant velocity. The filament that HH108MMS is forming within 
has a velocity gradient of 1.6 km s~^ pc~^ running from southwest to northeast. The CARMA 
N2H^ velocity map on the other hand shows significant structure on —10" scales within ±15" 
(—0.02 pc) of the protostar. Southeast of the protostar the emission is red-shifted and northwest 
there is blue-shifted emission, along the presumptive outflow axis and HCN emission mapped with 
CARMA (Tobin et al. 2011 in preparation). The outflow axis is determined from faint diffuse 
emission seen in Ks-band. The CARMA map also shows increased linewidth along the outflow; 
there is no indication of such an increase in the single-dish map, likely due to beam dilution. 
HH108MMS appears to be a prime example of the outflow impacting the envelope of a deeply 
embedded protostar. 



3.3.6. Serpens MMS3 
Serpens MMS3 is located within a complex network of filamentary structure in the Serpens B 



cluster ( jDjupvik et al.|[2006t [Harvey et al.||2006| ), shown in Figure [9j One of the most prominent 



filaments runs —0.1 pc in length into Serpens MMS3, see left panels of Figure [9) Directly west of 
the protostar the filament turns southward toward a small clustering of bright young stars. We 
also noticed that Serpens MMS3 has a faint companion separated by —7". 

The single-dish N2H^ map shows that the emission is highly pervasive throughout the region. 
The emission is peaked near the clustering of young stars in the southwest corner of the image 
in Figure [9| However, the emission is extended toward Serpens MMS3 and there is enhancement 
emission coincident with the large scale filament seen in 8/xm extinction. The VLA NH3 map 
reveals the structure of the region in substantially more detail. The interferometer resolved out 
whatever diffuse NH3 emission was in the region and the remaining emission directly correlates 
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to the highest extinction regions seen in the 8/xni image, with a peak coincident with the Serpens 
MMS3 protostar. The NH3 emission is stih extended toward the clustering of young stars in the 
southwest, but it is at the edge of the primary beam. 

The overah velocity structure of this region is confusing in the single-dish N2H^ velocity map 
because of the multitude of high-density structures in the region; however, there is a large scale 
gradient along the filament that Serpens MMS3 resides in and there is also a pocket of red-shifted 
emission next to the protostar. The VLA NH3 (1,1) map also shows this large scale gradient along 
the filament of Serpens MMS3 and red-shifted emission next to the protostar. The red-shifted 
emission appears over a region 15^' from the protostar and is extended in the direction of the 
outflow. We also detected an increased linewidth (^1 km s~^) next to the protostar. Directly 
southwest of the protostar the velocity gradient appears to resume the large-scale velocity trend 
exhibited northeast of the protostar. It is presently unclear if the kinematic structure near the 
protostar is related to the outflow or infall, but its proximity and extension along the outflow 
makes us suspicious. However, the broad linewidth is quite localized to the east of the protostar. 



3.3.7. HH211 



HH211MMS is a deeply embedded protostar on the outskirts of the IC348 cluster in the Perseus 



molecular cloud; emission from the central protostar itself only becomes evident at 70/im (Rebull 



et al.||2007 ) and has been found to be a proto-binary in the submillimeter ( Lee et ar||2009 ). We see 



a large absorbing structure in the 8/im extinction map shown in Figure 10, as well as its powerful 



outflow (McCaughrean et al. 1994; Gueth & Guilloteau 1999). The single-dish N2H+ emission 



associated with HH211MMS is very strong, peaked to the southwest of the protostar itself. The 
N2H^ emission also appears extended in the direction of the higher extinction areas. The N2H^ 
emission mapped with CARMA detects emission on small scales around the protostar, with the 
emission peak offset ^2" southwest of the protostar (not coincident with the single-dish N2H+ 
peaks) . The emission is more extended along the northwestern side of the outflow, consistent with 
the extinction seen in the 8/im image. 

The single-dish N2H+ velocity map shows a linear velocity gradient normal to the outflow and 
south of the protostar there is another velocity component in the dense gas. The transition between 
these two velocity components appears as an area of artificially large linewidth (an artifact from 
fitting); however, there are two sets of narrow emission lines present, not broad lines. The CARMA 
N2H^ velocity map also finds a linear gradient normal to the outflow as well as the second velocity 
component to the south. We also note that near the protostar the gradient is not perfectly linear 
at all scales. The deviance from a linear gradient is slight; however, it is present where we also have 



excellent signal-to-noise and this agrees with the velocity map by Tanner & Arce (2011). 



The linewidths in the single-dish data were quite low across the source, only 0.3 - 0.4 km s ^ 
with similar levels seen in the CARMA N2H+ map. We note that there is an area of increased 
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linewidth just southeast of the protostar, apparently at the base of the outflow. We suggest that 



the increased linewidth in this region is due to outflow interaction, in agreement with Tanner fc 



Arce (2011). In addition, the filament northeast of the protostar has a very narrow linewidth, <0.2 



km s ^, appearing both the single-dish and interferometer maps. 



3.3.8. IRAS 16253-2429 

IRAS 16253-2429 is a low- luminosity {LboI = O.25L0) Class protostar in the p Ophiuchus 
star forming region; it is also identified as Oph MMS 126 ( [Stanke et al.||2006 ). We noted in Paper 



I that this was one of the more "symmetric" envelopes seen in our 8/im extinction study. Its 



symmetric bipolar outflow has been traced in CO by Stanke et al. (2006) as well as in shocked H2 



emission from Spitzer IRS spectral mapping ( [Barsony et al.|[20T0 ). 



The single-dish N2H+ shown in Figure 1]_ correlates quite well with the 8/xm extinction. The 
emission peak is slightly offset from the location of the protostar and the N2H^ emission appears 
to be depressed at the location of the outflow cavities. The CARMA N2H+ emission shows similar 
features in that it strongly correlates with the regions of 8/xm extinction and there is less emission 
in regions occupied by the outflow cavities. The lack of emission within the outflow cavities is likely 
due to evacuation of envelope material and/or destruction of N2H+ by CO in the outflow (section 
4.2); there may also be some interferometric filtering-out of emission in this region. Furthermore, 
there appears to be a deficit of N2H+ emission near the protostar. 

The velocity field of the single-dish N2H+ map indicates that there is a very small velocity 
gradient across the envelope, approximately normal to the outflow; note the small velocity range 
occupied by the envelope and velocity gradient. The velocity map from the CARMA N2H+ data is 
more complex with several gradient reversals throughout the emitting region. However, the global 
gradient still seems to be present in interferometer data. Furthermore, a VLA NH3 map shows a 
velocity structure very similar to our CARMA N2H+ map (J. Wiseman, private communication). 

The single-dish N2H^ linewidth is quite small and constant across the envelope, whereas many 
other objects in our sample have linewidths which peak near the protostar. We also note that the 
linewidth is increasing in this source toward the edge of N2H^ emission along the outflow; this 
is likely an outflow interaction effect. The CARMA N2H+ map shows a similar small linewidth 
across the most of the envelope; however, there is a region of increased linewidth east and south 
of the protostar associated with an area of strong N2H+ emission; at this location there is a slight 
enhancement of linewidth in the single-dish map. 
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3.3.9. L1152 

The L1152 dark cloud is located in Cepheus, about 1.7 pc (20^) away from L1157 on the sky. 
LI 152 hosts three young stars; however, only one (IRAS 20353+6742) is classified as a Class 



object and it is the only one embedded in the main core of LI 152 (Chapman &; Mundy 2009) 



Paper I found that the main core of LI 152 appears to have a "dumbbell" morphology in which the 
northeastern core (see Figure [T2| ) appears to be starless and the southeastern core harbors IRAS 
20353+6742 (hereafter L1152). These two concentrations are connected by what appears to be a 
thinner filament of high density material. 



The single-dish N2H+ (J 



0) emission shown in Figure 12 exactly matches the morphol 



ogy of the extinction in the 8^m images. However, the peak N2H+ emission in the southwestern 
core is offset from the protostar by --20^' (6000 AU). The N2H+ map from CARMA shown in the 



right panel of Figure 12 only observed the southeastern core. The map confirms that the N2H+ is 
substantially offset from the protostar and there is no sub-peak at its location. However, the N2H+ 
emission appears to extend toward the protostar. 

The single-dish N2H+ velocity field exhibits a velocity gradient normal to the outflow of LI 152, 
noting that the protostar appears at the edge of the region exhibiting the gradient. The rest of 
the cloud, including the star-less core, has a fairly constant velocity; only varying by ^ 0.1 km/s. 
However, we do notice increased linewidths northeast and southwest of the protostar. Southwest 
of the protostar we can clearly see the jet from the protostar, possibly interacting with envelope 
material, then in the northeast there is nothing obvious happening at this linewidth peak in the 
8/xm image. However, the northeast linewidth peak is near the outflow axis and this could be the 
cause of the increased linewidth at this location. 

The velocity map from the CARMA N2H+ data tells a remarkably similar story to the single- 
dish data; the velocity gradient is only slightly better resolved. However, the most remarkable 
feature is in the linewidth map, where we clearly see an increase in linewidth on the axis of the jet 
that is visible in the S/xm maps. This appears to be a another very clear example of the outflow 
interacting with the envelope material, though the velocity field does not seem to show outflow 
effects. 



3.3.10. LI 527 



L1527 (IRAS 04368+2557) is an extensively studied protostar in Taurus. Benson &; Myers 



(1989) observed its compact NH3 core, from which Goodman et al. (1993) derived its velocity 



gradient. Subsequent observations indicated the possibility of infall in the envelope from H2CO 



observations by Myers et al. (1995). Furthermore, detailed modeling of its scattered light cavities 



observed in Spitzer IRAC imaging have been done by Tobin et al. (2008) and Gramajo et al.|p010 ). 
High-resolution mid-infrared imaging by Tobin et al. (2010a) found the signature of a large (R^200 
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AU) disk in scattered light. 

IRAC 8/xm imaging of this source revealed an asymmetric distribution of extinction, the north- 
ern side of the envelope is substantially more extended than the southern side. This asymmetry 
is also exhibited in our single-dish N2H^ shown in Figure [l3| in addition, the peak emission is 
also offset to the north of the protostar by -25'' (3500 AU). The VLA NH3 and CARMA N2H+ 
maps both show emission associated with the protostar, but the maps are somewhat difficult to 
interpret due to the likelihood of spatial filtering. In addition, the CARMA observation was done 
as a mosaic in order to cover the entire region of emission as the primary beam is only 70''; both 
maps seem to detect emission in the same general areas. 

The velocity field from the single-dish N2H^ map has a complicated morphology. There appear 
to be two velocity gradients in the map, one along the outflow (pointed out by Myers et al. ( 1995| )) 



and another normal to the outflow isolated by Goodman et al.| ( 1993| ). However, the gradient along 



the outflow is not linear, the velocities go from red to blue and back to red. The linewidth remains 
fairly constant throughout the map, with a minimum at the northeast and southwest edges of the 
map. 

The velocity fields from both the VLA and CARMA reveal further kinematic complexity in 
this system. We can see the consistency with the single-dish velocity map on large-scales; however, 
the N2H+ and NH3 maps show that there is a small-scale velocity gradient near the protostar. 
Notice that this small-scale velocity gradient is in the opposite direction as compared to the large- 
scale gradient. The linewidths of the N2H^ and NH3 exhibit a corresponding increase in the inner 
envelope, near these small-scale velocity gradients. This is the only protostellar envelope that where 
a velocity gradient reversal is seen going from large to small-scales. 



3.3.11. RNO43 

RN043 is protostar forming within the A Ori ring. On large scales the envelope is quite asym- 
metric, with several filamentary structures appearing to converge at the location of the protostar 



as shown in the left panels of Figure 14, RN043 also drives a powerful, parsec-scale outflow; CO 



emission has been mapped on small scales by Arce &; Sargent (2005) tracing an outflow cavity and 



on large-scales, tracing a ^5 pc long outflow (Bence et al. 1996). The N2H+ emission is mostly 



unresolved in the single-dish map as shown in Figure 14 The peak emission is located near the 



location of the protostar and there are slight extensions in the direction of the outflow. 

The CARMA N2H+ ( J = 1 ^ 0) map traces the small-scale structure seen in 8/im extinction 
very well. We also note a depression of N2H+ emission at the location of the protostar, consistent 
with observations of other protostars in our sample; see Section 4.5 for further discussion of this 
feature. East of the protostar there is a ridge of N2H^ emission which is composed of the three 
bright knots almost running north-south in the image extending ^35". The southern- most knot is 
associated with the highest column density region east of the protostar and the two northern knots 
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correlate well with an extinction filament running from the north into the envelope of RN043. 
In addition, this filament of 8/xm extinction and N2H+ emission are coincident with the brightest 
part of the outflow cavity in the Spitzer 3.6/xm image (Paper I). Directly west of the protostar, 
there is another peak of N2H+ emission and weaker N2H+ emission extended further west, in 
agreement with the 8^m extinction. Chen et al. (2007|) mapped this region in N2H+ using OVRO, 



the data agree quite well with our observations. However, our map appears to have recovered more 
large-scale emission, likely due to better uv-coverage at short spacings. 

The velocity field from the single-dish N2H+ map shows a large scale velocity gradient that 
is nearly normal to the outflow axis and there is an area of enhanced linewidth southeast of the 
protostar. The CARMA data reveal significant kinematic detail in the velocity field of the N2H+ gas. 
The CARMA velocity maps in Figure [M] clearly show red and blue-shifted sides of the envelope; 
however, separating those sides of the envelope is a sharp velocity jump from blue to red by 
^0.7 km s~^. Due to the overlapping lines at the location of the protostar, the N2H+ linewidth 
forms a line marking the jump in velocity. There also appears to be a north-south gradient in 



the interferometer data as well (the single-dish map hints at this). Chen et al. (2007) ignored 
the western, red-shifted portion of the envelope thinking that it was a line-of-sight alignment 
with another clump; however, the envelope has density increasing in Sjnm extinction toward the 
protostar on both sides (shown in Paper I), suggesting that the western side is indeed part of the 
same structure. 

We note that the most highly blue-shifted gas is not located directly adjacent to the protostar; 
this is likely due to the absence of N2H^ near the protostar, as mentioned earlier. Furthermore, 
small-scale emission of HCO^ was also detected with similar morphology to L1165 (Figure p|. 



Figure 15 shows that the centroid of blue and red-shifted emission are located normal to the outflow 
and are offset from each other by ^3^' (1400 AU). In RN043, the HCO^ line wings extend zb ^2 
km/s from the systemic velocity. If we assume that the HCO^ emission reflects only rotation, its 
velocities would imply an enclosed mass of 2.7 Mq. If only half of this velocity is due to rotation then 
the enclosed mass would be O.67M0. The bolometric luminosity of RN043 is ^S.OL©; comparable 
to LI 165 in both luminosity and mass. We have overlaid lines representing Keplerian rotation (or 



infall) on Figure 15, the O.67M0 curve matches the data much better that the 2.67M0 curve. 



Note that we have redefined the outflow position axis to be 20° east of north in contrast to the 
54° found by Arce &; Sargent ( |2005 ); our value is more accurate taking into account the outflow 



cavity observed by Spitzer (Figure 14) and CO maps from both [Arce fc Sargent (2005) and Bence 



et al. (1996). Furthermore, Chen et al. (2007) assumed the 54° outflow position axis, leading them 



to interpret the velocity gradient along the eastern ridge as symmetric rotation. The N2H+ gradient 
across the protostar has a very similar direction to the H^^CO^ and C^^O velocity gradients found 
by Arce &; Sargent (2005). However, our revised outflow axis and the observed N2H+ velocity 



structure, in conjunction with the H^^CO+ and C^^O data, lead us to suggest that we are likely 
not seeing envelope material being "pushed out" in this system, as suggested by |Arce fc Sargent 



(2005). Thus, the N2H+ velocity structure appears to reflect kinematic structure intrinsic to the 
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envelope. 



3.3.12. IRAS 04325+2402 

IRAS 04325+2402, sometimes referred to as L1535, harbors a multiple Class I protostellar 
system in the Taurus star forming region. The primary is possibly a sub-arcsecond binary with 
a wider companion separated by ^.2" (Hartmann et al. 1999[ ). The 8/xm extinction around IRAS 
04325 was found in our envelope study but not published in Paper I due to its low signal-to-noise; 



however, Scholz et al. (2010) noticed the 8/xm extinction in their study of the system. These authors 
pointed out that there is a bright diffuse region of emission at 3.6 and 4.5/im, at the peak of the 8/xm 
extinction. Furthermore, they noticed a dark band between the protostar and the 4.5/im diffuse 
emission, suggesting a dense cloud; however, there is a lack of 8/xm extinction at the location of the 
dark band. 

Our single-dish N2H+ map finds emission throughout the core surrounding the protostar with 
the peak emission coincident with the 8/xm extinction peak and the 4.5/xm diffuse scattered light 
peak. In fact, the N2H+ is peaked ^60^' northeast of the protostar, but the map does show a shght 
enhancement of N2H^ emission west of the protostar. Given this emission morphology, we suggest 
that the dark band see at 4.5/xm is really just a lack of material and that the diffuse emission is 
light from the protostar shining onto the neighboring star-less core. We have no interferometry 
data for this object, however we would not expect to observe substantial N2H+ emission peaked 
around the protostar, based on the single-dish map. 

The velocity structure of the N2H+ shows that there is a relatively smooth velocity gradient 
across the entire object with an increased gradient just southeast of the protostar. The linewidth 
of N2H+ however shows large increase along the outflow axis of the protostar. Thus, the outflow 
from the protostar may be interacting with the dense material in the surrounding core producing 
the increased linewidth. On the other hand, the velocity field does not seem to show effects from 
the outflow, similar to LI 152. 



3.3.13. 1483 



L483 is an isolated globule harboring a Class protostar ( [Tafalla et al.||2000| ). The envelope 
surrounding the protostar is quite large, ^0.15 pc in diameter with at least IO-2OM0 of material 
measured from 8/xm extinction in Paper I. The densest regions seen in 8/im extinction form a "tri- 
lobed" pattern that is also traced by 850/im emission ( j0rgense"ii||2OO4 ) . The single-dish N2H+ also 
follows this same pattern with the peak emission coincident with the protostar, see Figure pT} The 
VLA NH3 emission is not peaked on the protostar, but also follows the "tri-lobed" morphology. The 



N2H+ emission mapped with OVRO by J0rgensen (2004) is extended along the outflow; however 



this observation appears to have resolved out a significant amount of extended emission. 
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The velocity gradient from the single-dish N2H^ map is not normal to the outflow but is at an 
angle of ^45°. The VLA NH3 map shows a velocity gradient in the same direction as the single-dish 
data; however, the protostar is located in a pocket of blue-shifted emission. This is consistent with 



what J0rgensen (2004) observed. Furthermore, directly north of the protostar, there is an area 



of highly red-shifted emission seen in the NH3 map. The single-dish N2H+ also shows red-shifted 
emission in this region, but it is not as prominent due to the larger beamsize. This emission appears 
to come from another distinct velocity component in the cloud, as evidenced by the large linewidths 
in the NH3 map at the transition to the red-shifted emission. 

We also noticed that the N2H+ linewidth map shows a region of enhanced linewidth running 
across the envelope, nearly normal to the outflow. This region connects to where there is the second 
velocity component in the VLA NH3 map and this is also where the velocity field is most rapidly 
changing in the single-dish N2H+ map. Since the increased linewidth appears to be a global feature 
we do not attribute it to outflow effects and could be related to the initial formation of the dense 
core. 



3.3.14. L673 
The L673 dark cloud in the constellation Aquila has been the subject of a SCUBA survey by 



Visser et al. (2002) and two Spitzer studies by Tsitali et al. (2010) and Dunham et al. (2010). In 



Paper I, we highlighted a small region of the cloud exhibiting highly filamentary 8/xm extinction 



associated with L673-SMM2 as identified in Visser et al. (2002). There are more regions with 8/im 



extinction within the cloud that we did not focus on in Paper I, but are apparent in the images 



shown by Tsitali et al. (2010). The Spitzer IRAC data around L673-SMM2, show four point sources 



closely associated with the sub-millimeter emission peak and another 70/im source which may be 
a Class protostar ( jTsitah et aL]|2010D . 



The filamentary region around L673-SMM2 is shown in Figure 18, The N2H+ emission maps 
closely to the Sfim extinction and the N2H+ peak is centered on the small clustering of protostars. 
The peak NH3 emission from the VLA is located very near the N2H+ peak and the dense filament 
is further traced by the low-level NH3 emission; a substantial amount of extended emission is likely 
resolved out by the interferometer. 

The velocity field traced by the single-dish N2H+ appears to show a gradient along the filament 
going from north to south and there is an area of blue-shifted emission coincident with the southern 



protostar marked with an X in Figure [18) This southern- most protostellar source is comprised of 
three sources in higher-resolution Ks-band imaging (Tobin et al. in preparation). In the northeast 
part of the image, there is another velocity component of N2H+ present. The linewidths are fairly 
low across the filament with about a factor of two increase at the location of the protostars; there 
is an area of artificially large linewidth due to the second velocity component. 



The VLA NH3 map shows similar velocity structures that were present in the single-dish map; 
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however, it is now clear that the protostar near A5=0'' is located in an area of red-shifted emission 
while the southern protostar is still in a localized area of blue-shifted emission. The line-center 
velocity shift between these components is 0.4-0.5 km/s. The linewidth peak falls between the two 
main protostars, coincident with the region of peak NH3 emission. Also, the northernmost, deeply 
embedded protostar appears to be associated with a fairly ordered velocity gradient, north of the 
two more obvious protostars. 



3.3.15. L1521F 



L1521F is a dense core found in the Taurus star forming region. Bourke et al. (2006) found 



a deeply embedded protostar within what was previously considered a star-less core (Crapsi et al. 



2004, 2005). An approximately symmetric extinction envelope was found around L1521F, elongated 
normal to the outflow in Paper I. The N2H+ integrated intensity correlates very well with the 8/xm 
extinction. The NH3 observations from the VLA are also centrally peaked and show a flattened 
structure normal to the outflow. However, there is an extension to the east, along the outflow. 

The velocity structure of the core is complex and appears similar to that of LI 527. The 
N2H+ velocity field shows that there is emission blue-shifted relative to the protostar normal to the 
outflow. Along the outflow there is red-shifted emission toward the edge of the envelope. [Crapsi 



et al. (2004) examined the velocity structure of L1521F finding that the average gradient across 
the core was 0.37 km s~^ pc~^ with a position angle of 180°. With our two dimensional fitting 
we derive a gradient of 0.76 km s~^ pc~^ and a position axis of 239°. The differences between 
out results likely come from mapping a larger area around the core, which detects more red-shifted 
emission in the western side of the map, influencing the gradient fit. Otherwise, the emission and 
velocity structure are quite similar. 

The NH3 velocity map shows similar structure to the N2H+ map, and near the protostar 
there appears to be a gradient emerging normal to the outflow on small-scales. However, the 
NH3 emission is optically thick toward the center of L1521F: the satellite and main lines have 
approximately equivalent intensities. Therefore, we cannot obtain a better measure of the small- 
scale kinematic structure. The N2H+ linewidth map shows a roughly constant 0.2 - 0.3 km/s 
linewidth across the map. The NH3 linewidth is similarly low, except at the southern end of the 
envelope where the blue-shifted emission is present. 



3.3.16. Perseus 5 
Perseus 5 is a relatively isolated core in the Perseus molecular cloud, just northeast of NGC1333 



and observed by Caselli et al. (2002). The protostar is deeply embedded and is obscured shortward 
of S/xm with only its outflow as a prominent signpost. It was discovered to have an asymmetric 
extinction envelope around it in Paper I. We did not have the opportunity to take single-dish N2H+ 
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observations of this object, but we did take data with CARMA as shown in Figure 20, The N2H+ 
intensity image shows that the entire extinction region is not well traced by the interferometric 
N2H+. The data indicate that substantial emission around this source is resolved-out indicated 
by the strong negative bowls in the image. However, NH2D (another molecule we observed) does 
seem to fully trace the envelope seen in 8^m extinction, since emission from this molecule is more 
spatially compact and not filtered-out by the interferometer. 

The velocity structure is complex in both N2H+ and NH2D, showing a blue-shifted feature east 



of the protostar along the outflow (Figure 20). Furthermore, both tracers show a similar gradient 
along the outflow direction; NH2D shows increased linewidth through the envelope, close to the 
outflow direction and there are several regions of enhanced linewidth in N2H+ along the outflow. 
Furthermore, there may be a gradient normal to the outflow as seen in both the N2H+ and NH2D 
velocity maps. However, the outflow seems to be significantly influencing the kinematics of the 
dense gas. 



3.3.17. IRAS 03282+3035 
IRAS 03282+3035 is an isolated, deeply embedded Class protostar located in the Bl-ridge of 



the Perseus star forming region (J0rgensen et al. 2006). Mid-infrared emission from the protostar 



itself is quite faint, but appears as a point-source at 8/xm and Chen et al. ( 2007] ) identified it as a 
binary in millimeter continuum emission. The IRAC 8/xm extinction toward this object in Paper I 
highlights a rather complex morphology on large scales; however, near the protostar the extinction 
appears to be concentrated into a filamentary structure. 



The single-dish N2H^ observations in Figure 21 trace the large-scale extinction morphology 
very well and the emission is observed to be quite extended, with the emission peak slightly offset 
from the protostar along the outflow. The N2H+ emission also ends at the northeast edge of the 
core where the extinction rapidly falls off. The VLA NH3 map traces a filamentary structure on 
large-scales north and south of the protostar. Furthermore, the emission double-peaked, with the 
individual peaks located north and south of the protostar, in agreement with the N2H+ emission 



shown by Chen et al. (2007). 



The velocity field derived from the single-dish N2H^ map shows a strong velocity gradient in 
the direction of the outflow; however, there also appears to be another gradient that is normal to 
the outflow on the southeast side of the envelope. We also note that there is a strong linewidth 
gradient in the direction of the outflow (same direction as the line-center velocity gradient) with 
the largest linewidths appearing on the west side of the envelope. The NH3 velocity map from 
the VLA again finds the velocity gradient in the direction of the outflow along with blue-shifted 
emission north and south of the protostar. The linewidth of the NH3 emission is peaked just north 
and south of the protostar indicative of dynamic motion in the line of sight. 



Chen et al. ( |2007[ ) suggested that the outflow may be interacting with the envelope, causing 
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the gradient along the outflow. However, the IRAC data show that the outflow is well-collimated 
and appears only able to affect the material extended to the northwest of the protostar. To further 
examine the kinematic structure of this system, we have overlaid the NH3 channel map contours on 



the CO map from Arce & Sargent (2006) in Figure 22, The top and bottom rows show the portions 



of the NH3 lines that are not blended while in the middle panel the two lines are blended. At 
red-shifted velocities in the top panel, notice how the different parts of the emission come into view 
going toward lower velocities. The top right panel shows the red-shifted emission being emitted 
from two very thin structures running north-south. In the bottom panels of Figure 22, we can 
see that the more blue-shifted emission is still confined to a very thin structure running north- 
south; however, its location has shifted further to the east compared to the upper panels. Thus, 
we suggest that rather than an outflow interaction giving rise to the global velocity structure, a 
more probable explanation is that the envelope is filamentary and that the velocity shifts are flows 
along the filament toward the protostar. However, certain regions of the envelope that are spatially 
coincident with the outflow cavities may be affected by the outflow. 



3.3.18. HH270 VLAl 

HH270 VLA 1 is located in the L1617 cloud near Orion and is the driving source of the HH270 
outflow (Reipurth et al. 1996; Rodriguez et al.|[l998D . This outflow in particular appears to be 



deflected and colliding with a neighboring core (the driving source of HHllO (Rodriguez et al 



19981)). In Paper I, we found that this protostar exhibits 8/xm extinction in a filamentary envelope 



that is extended along the outflow and only on one side of the protostar. Furthermore, it was noted 
that the 4.5/xm scattered light is brighter at the edge of the extincting envelope, indicating that 
extinction structure is indeed extended along the outflow and not due to a complex line of sight 
projection effect. The single-dish N2H+ emission from this source is peaked on the 8/xm extinction 
with no peak coincident with the protostar. The CARMA N2H+ map is still tracing the 8/im 
extinction and shows a sharp decline of emission at the location of the protostar; the emission still 
appears one-sided even at high resolution. 

The single-dish N2H+ velocity map shows a gradient along the outflow with a linewidth en- 
hancement just southwest of the protostar. This region of high linewidth is likely an unresolved 
velocity gradient, which is then resolved in the CARMA velocity map. The pattern in the CARMA 
velocity map closely matches the shape of the scattered light outflow cavity visible at 8^m in Figure 
[23} We suggest that these features are due to an outflow interaction with the surrounding material. 
We note that the data do not show a linewidth increase; therefore the outflow seems to be inducing 
bulk motion, but not an increase in linewidth. 
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Discussion 



4.1. Large-Scale Velocity Field 

The large-scale velocity gradients in dense molecular cores have long been interpreted as ro- 



tation (possibly solid-body) (e.g. Goodman et al. 1993; Caselli et al. 2002; Belloche et al. 2002 



Belloche fc Andril|2004[ Chen et aL||2007 ). If the large-scale motions are indeed rotation, then the 
circumstellar disk forming from the collapse of the surrounding envelope should be rotating in the 



same direction as the surrounding cloud, assuming axisymmetric collapse (e.g. Shu et aL]|1987 ; Bo- 
denheimer||1995 ) . However, if the core has fragmented, then the angular momenta of the individual 
disks could be misaligned, but the sum of the angular momenta would still be aligned with the 
cloud rotation. The direction of the jet and outflow of the system is commonly thought to reflect 
the current angular momentum vector of the protostar and disk system. This is because outflows 
and jets are thought to be related to the magnetic fleld of the protostar and disk which should 



be along the rotation axis (Pudritz &; Norman 1983; Shu et al. 1994) and resolved observations 



of disks have shown that jets are oriented normal to the disk midplane. Thus, we have used the 
outflow direction and protostar locations as guideposts as to where we should measure the velocity 
gradients relative to and in which direction because the equatorial plane of the envelope is where 
we are most likely to observe rotation. 

The 2D fltting method (see Section 2.6.2) shows that eleven out of sixteen systems have a 
gradient direction that is within 45° of normal to the outflow axis in the single-dish data; however, 
there is considerable spread in this distribution. We caution that the PA of the gradients calculated 
by this method may have systematic error due to the velocity flelds not being uniform and the 
velocities are sampled over a region that is not symmetric, leading to a bias of data points in a 
particular direction. Visual inspection of the velocity flelds does yield a similar result to the 2D 
fltting, with eleven normal to the outflow and flve not normal; IRAS 03282+3035, HH270 VLAl, 
L1527, L1521F, and L483 were not normal to the outflow. In contrast to the 2D fltting, we visually 
identify L483 and HH270 VLAl as not normal to the outflow and LI 157 and Serpens MMS3 to be 
normal. We further caution that most velocity flelds have structure not easily described by a single 
position angle or velocity gradient. Nevertheless, as a simple signiflcance test of the gradient-outflow 
direction relation, we can apply the binomial distribution if we consider the <45° and >45° as two 
bins and that the gradient direction relative to the outflow may be oriented randomly between 



and 90° with a mean of 45° (Bevington 1969). Thus, the probability of the relative directions 



falling within <45° or >45° should be equal. The chance of eleven or more objects out of sixteen 
with randomly oriented velocity gradients falling within the >45° bin is only ^10%. 

We also attempted to see if there is any trend with envelope mass and velocity gradient, 
we compared the observed velocity gradients and linewidths with the escape velocity in Figure 
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The envelope mass is calculated from the dust mass measured from Sfim extinction (Paper 
I, Table 1) plus 1 Mq to account for the central object; note that the central object masses 
may be overestimated while the envelope masses are likely underestimated (Paper I). There is 



29 



no apparent trend of velocity gradients or linewidth with increasing envelope mass (in terms of 
escape velocity). This figure further demonstrates that these protostellar systems are consistent 
with being gravitationally bound and are not supported by rotation or turbulence, consistent with 
many previous studies (e.g. Goodman et al.|[l993 ; Caselli et al.||2002 ; Chen et "aL||2007 ). 



Furthermore, the smallest velocity gradients in our sample appear to be coming from the most 
symmetric envelopes L1521F, IRAS 16253-2429, and L1157, and the largest gradients are found 
in the morphologically complex HH211 and IRAS 03282+3035 systems. This could be taken to 
mean that the velocity structure is no strongly influenced by complex projection effects in the 
symmetric systems and are observing slow envelope rotation, while the more complex systems have 
projection effects altering the observed velocity structure. However, despite this trend, we note 
that the envelope of the low-luminosity source IRAM 04191 is approximately symmetric (Paper I) 



and has a velocity gradient of 17 km s ^ pc ^ (Belloche et al. 2002; Belloche &; Andre 2004), a 
clear counter-example to this morphological trend. 

The simple interpretation of our velocity gradient data as a whole is that we are observing core 
rotation in these systems and that they have a variety of angular momenta ; some angular momenta 



being quite large (Figure 25). However, if we were observing pure rotation on large-scales, then one 
would expect the velocity gradient directions to be much more clustered toward being orthogonal 
to the outflow rather than the broad distribution shown in Figure [27| Therefore, either the complex 
morphology makes the rotation ambiguous, causing the velocity gradients to not be orthogonal to 
the outflow or we are not observing pure rotation. In either case, the angular momenta derived 
from the velocity gradients will be suspect at best. Given these complicating factors, we regard the 
velocity gradients taken as ID cuts normal to the outflow as most likely to be probing the velocity 
gradients due to rotation, but these values should be regarded as upper limits. 

If we are not observing pure rotation, the only other dynamical process which should give rise 
to ordered velocity structures are infall and/or outflow entrainment. The outflows from the proto- 
stars in our sample are highly collimated, limiting their ability to affect the large-scale kinematic 
structure, see Section 4.3 for further discussion, which leaves infall as the only other mechanism to 
contribute to the velocity field. This would require that infall is happening on large-scales and that 
collapse would need to be outside-in and not inside-out. Large-scale infall is shown to be possible 



in numerical simulations of complex, filamentary cores forming within a molecular cloud by Smith 



et al. (2011). Their simulations show that there is infall onto filaments (scales of 0.1 pc to 0.01 pc) 
from the surrounding molecular cloud, with subsequent infall from the filament to the protostar 
(sink particle) (< 0.01 pc). Thus, infall convolved with rotation in real envelopes could lead to 
some of the complex velocity structures that we observe. 

In summary, it is difficult to interpret the large-scale velocity gradients wholly being due to 
rotation since the envelopes are asymmetric and the gradients are often not normal to the outflow. 
Furthermore, flows or infall along a filament could give a similar signature to rotation simply from 
geometric projection, and if an envelope is filamentary, infall along the envelope would also produce 
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a velocity gradient normal to the outflow. We therefore suggest that a component of infall velocity 
could be projected along our line of sight, entangled with the rotation velocity, resulting in the 
large velocity gradients. This issue of rotation versus projected infall will be further explored in an 
upcoming paper. 



4.2. Small-scale Velocity Structure 

The high-resolution interferometer data are essential for probing the kinematics of the envelope 
at scales smaller than 5000 AU. These data enable us to localize the gas in the envelope and to assess 
the dynamical processes at work. Small-scale velocity sub-structure, beyond an ordered/linear 
velocity gradient, is found in the interferometric observations of most envelopes, including LI 157, 
L1165, HH108IRS, Serpens MMS3, L1527, CB230, IRAS 03282+3035, and RN043. These features 



are sometimes apparent in the velocity maps or profiles shown in Figure [24| but they also appear 
as increased linewidth in the inner envelope. The small-scale velocity structure is generally found 
on ^2000 AU scales in the envelopes. This radius could possibly be the centrifugal radius where 
material can be rotationally supported against gravity and the increased rotation velocity makes 
this region stand out against the rest of the envelope in velocity. However, assuming the large-scale 
velocity gradients reflect rotation, there would not be enough angular momentum for material to 
be rotationally supported at this radius. 

We found that the relationship between outflow axis and velocity gradient direction in the 
interferometer data trends even more strongly toward being normal to the outflow that in the 
single-dish data. Twelve of fourteen systems have gradients within 45° of the outflow in the right 



panel Figure [27| as found by the 2D fitting technique. Using the same statistical analysis as the 
single-dish gradients, this result being due to chance is ^0.6%. Note that this plot is missing L1157 
and Serpens MMS3 since their complex velocity fields could not be reliably fit. Visually, we find 
that twelve envelopes clearly have gradients normal to the outflow and four do not (HH108MMS, 
L483, IRAS03282, and HH270VLA1); L483 and HH270 VLAl were counted as within 45° of normal 
to the outflow. Note that RN043, Perseus 5, and L1521F appear to have gradients normal to the 
outflow, but they also have complex velocity fields. L1157 is left out of this analysis because the 
gradient directions from the CARMA and VLA data differ by ^80°. 

The velocity gradient direction is generally consistent between both large and small scales in 



the single-dish and interferometer data as shown in Figure 28 This could mean that the dynamical 
processes observed at large-scales are also responsible for the kinematics observed at small-scales. 
The protostellar systems which show substantial deviation from large to small-scales are LI 527 and 
L1521F. The difference in L1527 marks a velocity gradient reversal from large to small-scale (as 



also shown in Figure 24) and the difference in L1521F reflects that the kinematic structure was 
not well-resolved in the single-dish data. The N2H+ and NH3 lines in L1521F are also extremely 
optically thick making the velocity field derived from fitting uncertain since we cannot probe all 
the gas along the line of sight. 
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The preference for the vast majority of velocity gradients to be within 45° of normal to the 



outflow contrasts with Chen et al. (2007), where only two of nine targets had this feature. Moreover, 
in Volgenau et al. (2006), only one protostellar core of three observed had a well-ordered velocity 
structure. This result may be due to environment since the ordered structure was found in L1448 
IRS3, a more isolated system, and the other sources were located in the more complex environment 



of NGC1333. Our sources in contrast are generally isolated, like those in Chen et al. (2007). As 
noted in the preceding paragraphs, we do see several objects with similarly complex velocity fields 
in our data; however, our greater number of observed systems likely enabled us to find more with 
ordered velocity fields. 

The velocity gradients fit for the interferometer data are systematically larger than those found 
for the single-dish data, with the average being 8.6 km s~^ pc~^ from ID fitting, the distribution 
is shown in Figure [26} This factor is nearly equal to our increase in resolution in the interferometer 
data as compared to the IRAM 30m, but the object-to-object increase is more varied. Some of the 
small-scale velocity structure is due to outfiow interactions in the inner envelope as we will discuss 
in more detail in the following section; however, we have attempted to mask out these regions when 
fitting the velocity gradients. Like the single-dish data, we regard the ID fitting method to more 
accurately refiect the velocity gradient intrinsic to the envelope itself since the 2D method will be 
more susceptible to outfiow effects on the envelope kinematics. The only objects not showing large 
increase in the velocity gradient on small-scales are HH211 and IRAS 16253-2429 whose small-scale 
velocity gradient may be slightly overestimated. 

We strongly cautioned in the previous section about interpreting the velocity gradients as ro- 
tation on large-scales and we are again hesitant to interpret the small-scale gradients as rotation 
and/or spin- up. This is because most envelopes are highly filamentary and on small-scales pro- 
jections effects on the velocity structure will be even more apparent since both infall and rotation 
velocities increase at small radii. The velocity fields themselves on small-scales are not well-ordered 
as one might expect from rotationally dominated motion. Highlighting a few examples from Section 
3.2: Serpens MMS3 shows deep red-shifted emission and increased linewidth only on one side of the 
protostar, the velocity gradient in HH108IRS reverses itself just past the protostar, the gradient in 
L1527 on small-scales is opposite of large scales, and RN043 has an abrupt velocity jump across 
the envelope. The complex nature in the velocity fields of many sources do not necessitate an 
interpretation as rotation. Furthermore, infall velocities will always be necessarily be larger than 
rotation since the envelopes are not rotationally supported ( Chen et"aL]|2007 ). Thus, the only way 



to robustly separate infall velocities from rotation is at < 1000 AU scales, near the centrifugal 
radius where material can be rotationally supported. 

Two systems in our sample (LI 165 and RN043) show high- velocity line wings in the inner 
envelope that could indeed refiect significant rotation. In both sources, the emission appears to 
be coming from a radius of ^600 AU. Our data indicate that the observed velocities could refiect 
rotationally supported motion around ^O.5M0 central objects. Higher resolution and higher signal- 
to-noise data are needed to accurately centroid the high velocity emission in order to more precisely 
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constrain the enclosed masses. In the future, ALMA could be used to spatially resolve the high- 
velocity emission in order to fully separate rotation from infall. Further analysis and modeling 
of the velocity structures observed with the interferometric data with the goal of determining the 
kinematic processes at work in the envelopes will be presented in an upcoming paper (Tobin et al. 
2011 in preparation). 



4.3. Outflow Induced Kinematic Structure 



Our knowledge of protostellar outflows has been greatly enhanced in recent years (e.g. Bachiller 



1996; Arce et al.|2007, and references therein) and their possible effects on the surrounding envelope 



have been characterized (Arce fc Sargent|2006|). Furthermore, IRAC imaging from the Spitzer Space 



Telescope, in addition to near-IR imaging from the ground, can give a strong constraint on the 
outflow axis and cavity width (e.g. Scale fc Looney|[2008 ) . These data complement observations of 



outflow tracers such as CO and give a more complete picture of how the outflow may be impacting 
the protostellar envelope. 

Several protostars in the sample have velocity structures that are strongly suggestive of outflow 
effects on the N2H+ and NH3 gas kinematics. L1157, HH108MMS, and Perseus 5 appear to be 
significantly affected by the outflow, while HH108IRS, HH211, IRAS 03282+3035, L1152, and IRAS 
04325+2402 only appear to be mildly affected. In the mildly affected cases, the large-scale bulk 
motion of the envelopes (line-center velocity) does not appear to be affected by the outflow, rather 
we see the effects in the linewidth at both large-scales and small-scales. This probably means that 
only a small portion of the total cloud mass is being affected by the outflow. The bulk motion 
effects are revealed in the strongly affected cases, generally on small-scales and only visible in the 
interferometer data. 



Arce &; Sargent (2006) presented an empirical model for how the outflow will affect the envelope 
during protostellar evolution, concluding that the outflow is ultimately responsible for disrupting 
the protostellar envelopes. The number of objects we find showing outflow effects in the kinematic 



data strongly support the outflow-envelope interaction framework put forward by |Arce fc Sargent 



(2006) and enable us to offer some further input to this empirical framework. 



Perseus 5 and HH108MMS appear to be some of the youngest objects in our sample, as evi- 
denced by their deeply embedded nature and lack of visible outflow cavities in 3.6/xm or Ks-band 
imaging. In these systems, the outflow seems to be having the greatest impact on the kinematics 
in both line-center velocity and linewidth. Therefore, the effects of the outflow on the kinematic 
structure may be most prominent during its initial breakout of the envelope, early in the Class 
phase. Smaller-scale effects on the envelope can clearly be seen in the case of LI 157 where the 
envelope material may be entrained at 1000 AU scales. Thus, the outflows could be carrying sig- 
nificant momentum at wide angles near the protostar in order to be actively forcing inner envelope 
material out. If the outflows are impacting the envelope at small-scales and wide angles, an impor- 
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tant question remains as to how much the outflow can quench infall in filamentary envelopes. In 
L1157, the amount of entrained gas appears tenuous, only ^0.04 Mq out of 0.7 Mq in the inner 
envelope, following a simple analysis ( Goldsmith fc Langer||l999 ) assuming an N2H+ abundance of 



10~^ per H2. On larger-scales (> lOOOAU), we expect that material extended normal to the outflow 
in filamentary envelopes will not be strongly influenced by the outflow, as the outflows appear to 
stay well collimated throughout the Class phase. 



4.4. Linewidths 

The linewidths in the envelopes are generally quite small away from the protostar. In most 
cases, the linewidths are 0.2 - 0.5 km s~^ in the single-dish data; linewidths averaged over the 
entire source are given Table 3. Note that these linewidths are substantially broader than the lOK 
thermal linewidth of N2H+ which is ^0.13 km/s, and is generally attributed to turbulent motions 
and/or unresolved velocity gradients along the line of sight. However, these linewidths are not large 



enough to make the envelope unbounded, as shown Figure 29 



The interferometer data find even smaller linewidths toward the outer edges of the envelopes; 
there were several cases where the lines are less than 2 channels wide. The narrow line-widths in 
the interferometer data likely stem from the larger scale structure being resolved out; the larger 
scale emission may have some turbulent or large-scale infall velocity component. The regions of 
broad linewidth in the interferometer observations appear to be directly related to the outflow, 
increased line-of-sight motion, and/or heating from the protostar; not a transition to a turbulent 
core (cf. |Chenet aLl[2007| ). 



NH3 linewidths were used by Pineda et al. (2010) to probe the kinematics of the large-scale 



molecular core, detecting a transition from the quiescent core to the turbulent cloud. We looked 
for such an effect in our sensitive N2H+ data but did not find similar structure. We also did not 



detect N2H+ emission on the scales for which Pineda et al. (2010) were able to detect NH3, despite 



our high sensitivity. We believe that this is likely due to the differing critical densities between 
N2H+ and NH3 (-2x10^ cm-^ for NH3 (1,1) and -1.4x10^ cm-^ for N2H+ ( J = 1 ^ 0)). 

We did, however, find that some maps (i.e. L673, HH211, HH108, RN043) had regions with two 
distinct velocity components. The regions of overlap between the components appear as artificially 
large linewidths in the maps generated by the hyperfine fitting routine (section 2.6.1). We show 
three N2H+ spectra from HH211 taken at three positions showing the different velocity components 



in Figure [30t this is indicative of what takes place in the other regions showing this feature. Notably, 
the additional velocity components tend to appear toward the edge of the maps, except in RN043 
where the transition takes place near the protostar. The distinct velocity components are located 
in regions that appear form a contiguous structure when viewed in 8/xm extinction. The second 
velocity component generally appears about ^0.05 pc from the nearest protostar. Thus, the reason 
for multiple components could have to do with the initial conditions of the clouds themselves. 
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Recent simulations have suggested that cohiding clouds could be an important component to setting 



up the initial conditions for star formation (e.g. Heitsch et al. 2006). 



4.5. Chemical Effects on Molecular Tracers 

The kinematic data presented are based on the molecular tracers N2H+, NH3, and HCO^. We 
have used N2H+ and NH3 relatively interchangeably, since they appear to trace the same kinematics 



and physical conditions (see Section 3.1 and [Johnstone et aL] ( |2010 )). This makes sense because we 
know that in the pre-stellar phase, the formation of N2H^ and NH3 appear to be linked given their 
similar abundance distributions and only deplete onto dust grains at very high densities ( [Bergin fc 
TafallapOOT ) . In many of our observations, N2H+ and NH3 only appear to trace the gas on scales 
>1000 AU from the protostar; the emission generally peaks near the protostar, but not directly 
on it. This indicates a drop in abundance either due to depletion or destruction of the molecules. 
Several of the protostars for which there are both interferometric NH3 and N2H+ observations (in 
this paper or in the literature) show decreased emission at the location of the protostar in both 



N2H+ and NH3. Observations of CB230 and IRAS 03282+3035 in NH3 (Figures ^ and 21) and 



N2H+ (Chen et al. 2007) are clear examples of NH3 and N2H+ not being peaked coincident with 



the protostar. LI 157 also exhibits this effect in N2H+ from the PdBI observations and also does 
in NH3 if the longer baselines are given more weight in the VLA map. Finally, the low-luminosity 



source IRAM 04191 shows a similar depletion pattern in both of these tracers (Belloche &; Andre 



2004, J. Mangum, Private Communication). 



We can understand the decrease in N2H+ emission in terms of molecular destruction by re- 
actions with other molecules. In the inner envelope, where the temperatures rise above 20K, the 
CO that has depleted onto the dust grains ( Bergin et al.||2002 ) is released back into the gas phase. 
CO and N2H+ rapidly react to form HCO^; this is the dominant destruction mechanism for N2H+ 
( Aikawa et al.||2001 ; Lee et al.]|2004| ). This is the same reason that N2H+ is not seen in outflows and 
why the N2H^ is often not centrally peaked on the protostars in our interferometric observations. 
NH3 however, is not directly destroyed by CO, but rather HCO+ ( Lee et al.|2004 ). HCO+ will read- 
ily form within the region of CO evaporation making it available to react with NH3. Alternatively, 
NH3 could also become depleted onto dust grains in an ice mantle, and if it is well-mixed with the 
water ice, then it would only be evaporated at temperatures >100K. NH3 ice is frequently observed 
in the envelopes surrounding protostars ( Bottinelli et al.|[2010 ) via mid-infrared spectroscopy. The 
absorbing NH3 ice should be in the inner envelope since NH3 only depletes onto grains at high 
densities. If any of the NII3 ice were released into the gas phase, then HCO+ would be present to 
destroy it. 

Since both NH3 and N2H+ are not present on scales <1000 AU, we must look for other tracers 
to probe the kinematic structure on these scales. In two protostars, we have been able to use HCO^ 
to trace the small-scale kinematic structure. The HCO^ emission from LI 165 and RN043 show 
high-velocity wings on small-scales inside the innermost N2H+ emission; the centroids of the red 
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and blue-shifted emission are offset from the protostar, normal to the outflow direction, at radii of 
-600AU. 



Chemical models of an infalling protostellar envelope have been calculated by Lee et al. (2004), 
showing that HCO^ becomes enhanced at small-scales after formation of the protostar. This reflects 
the evaporation of CO ice, releasing CO back into the gas phase where HCO^ is readily formed. In 
addition, the primary destruction pathway of N2H^ from CO results in the formation of HCO^ and 
N2. Thus, we can see why the N2H^ is tracing the gas at larger radii and lower velocities, while in 
the inner envelope HCO^ is readily being formed and can then trace the small-scale high-velocity 
gas. 

Our dataset explicitly shows how multiple tracers can be used to gain a more complete picture 
of the kinematics in protostellar envelopes. N2H+ and NH3 are excellent tracers of the cold, dense 
gas on scales from -1000-2000 AU out to -10000-20000 AU. Inside of lOOOAU, an abundant tracer 
of the warm, inner envelope tracer must be observed; HCO^ works quite well in two out of nine 
sources observed with CARMA and other tracers and/or higher-J transitions of HCO^ may work 
as well (see Lee et al. (2009) and Brinch et al. ( 2007[ )). However, these other tracers are often found 



in outflows and they must be observed with sufficiently high resolution to confirm their origin in 
the envelope and not the outflow. In the future, ALMA may be able to observe inner envelope 
tracers to resolve the motion of the dense gas in the inner envelope, tracing infall onto the disk. 



5. Summary 

We have conducted a single-dish and interferometric survey mapping emission of the dense gas 
tracers N2H+ (J = 1 ^ 0) and/or NH3 (1,1) in envelopes around low-mass protostars. Many of 
these envelopes are known to be morphologically complex from 8/im extinction mapping. We used 
these data to map the line-center velocity and linewidth across these envelopes. We quantitatively 
measured the velocity gradients and their directions in order to characterize the dense gas kinematics 
in the complex protostellar envelopes, our specific results are as follows. 

1. Ordered velocity fields are present on large-scales in most protostellar envelopes from the 
single-dish sample. In eleven out of sixteen cases, the velocity gradients appear to be within 45° of 
normal to the outflow axis with an average gradient of ^ 2.3 km s~^ pc~^ , depending on fitting 
method. The velocity gradients could be due to core rotation; however, the velocity gradient 
position angles do have a broad distribution with respect to the outflow direction. Furthermore, 
the strongly asymmetric nature of the envelopes and the fact that most envelopes in the sample 
are substantially bound leads us to suggest that we may be seeing a component of infall projected 
along our line of sight entangled with rotation. We find evidence of multiple components of N2H+ 
emission in several clouds, possibly relating our observations to a colliding-cloud formation scenario. 

2. The small-scale kinematic structure observed by the interferometers appears to be gravita- 
tionally dominated by the central protostar and it likely originates from a combination of infall and 
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rotation in some cases. The average velocity gradient in the interferometer data is 8.6 km s~^ pc~^ , 

o 

with gradient directions within 45 of normal to the outflow in twelve out of fourteen cases. The 
complex velocity fields in many systems suggest that interpreting the velocity gradients as pure 
rotation is incorrect. Only on the smallest scales (<1000 AU) will rotation differentiate itself from 
infall. Furthermore, multiple tracers must be used to gain a complete picture of the kinematic 
structure of the envelope down to sub- 1000 AU scales due to depletion of the cold gas tracers 
within -1000 AU. In the cases of RN043 and L1165, we were able to use HCO+ ( J = 1 ^ 0) to 
trace inner envelope kinematics on scales of —600 AU. 

3. Outflows do impact the envelope kinematics derived from NH3 and N2H+ in some systems 
and their effects are most prominent at —1000-2000 AU scales. Clear outflow effects on the large- 
scale kinematics in the single-dish data are only seen for five systems, but the only evident kinematic 
effect is increased linewidth along the outflow. The effects are most pronounced in the interferometer 
data for the deeply embedded sources HH108MMS and Perseus 5. Furthermore, we see a remarkable 
case of the outflow possibly entraining the inner envelope of LI 157 on —1000 AU scales. 
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A. Continuum Data 

In conjunction with our CARMA and PdBI N2H+ observations, the 3mm continuum emission 
was observed for those 11 sources. The continuum data are overlaid on the IRAC 8/im images 



for ah but HH108IRS/MMS (we shown 70/xm data due to lack of IRAC imaging) in Figure 31 



The 3mm continuum sources are point-like in most cases, HH108IRS and HH211 appear to show 



37 



extended structure at the 3cr level, consistent with envelope dust emission. LI 157, LI 165, and 
HH270 are also slightly extended, but in the direction of their outflows, indicating that there may 
be a component of free-free jet emission or heated dust along the outflow in the 3mm continuum 
data. In all cases, the continuum sources are coincident with the 8/im/70/xm point sources. There 
are a few cases (L1157, HH211, and RN043) where there is not a clear point source at S/xm and 
the 3mm emission is located between the outflow cavities where there is significant extinction. In 
these cases, the 24/im point source is then found to peak at the same location of the 3mm source. 
We note that the protostar IRAS 16253-2429 shows a shght offset between the S/xm point source 
and the continuum detection; however, the detection is barely 3cr, thus we do not believe that this 
offset is real. The continuum fluxes at 3mm are listed in Table 11. 



REFERENCES 

Aikawa, Y., Ohashi, N., Inutsuka, S., Herbst, E., & Takakuwa, S. 2001, ApJ, 552, 639 

Andre, P., Ward- Thompson, D., & Barsony, M. 1993, ApJ, 406, 122 

Arce, H. G., & Sargent, A. I. 2005, ApJ, 624, 232 

— . 2006, ApJ, 646, 1070 

Arce, H. G., Shepherd, D., Gueth, F., Lee, C., Bachiller, R., Rosen, A., & Beuther, H. 2007, Protostars and 
Planets V, 245 

Arquilla, R., & Goldsmith, P. F. 1986, ApJ, 303, 356 

Bachiller, R. 1996, ARA&A, 34, 111 

Barsony, M., Wolf-Chase, G. A., Ciardi, D. R., & O'Linger, J. 2010, ApJ, 720, 64 

Belloche, A., & Andre, P. 2004, A&A, 419, L35 

Belloche, A., Andre, P., Despois, D., & Blinder, S. 2002, A&A, 393, 927 

Bence, S. J., Richer, J. S., & Padman, R. 1996, MNRAS, 279, 866 

Benson, P. J., & Myers, P. C. 1989, ApJS, 71, 89 

Bergin, E. A., Alves, J., Huard, T., & Lada, C. J. 2002, ApJ, 570, LlOl 

Bergin, E. A., & Tafalla, M. 2007, ArXiv e-prints, 705 

Bevington, P. R. 1969, Data reduction and error analysis for the physical sciences, ed. Bevington, P. R. 

Bodenheimer, P. 1995, ARA&A, 33, 199 

Bonneh, I. A., & Bate, M. R. 1994, MNRAS, 269, L45 

Bottinelli, S., Adwin Boogert, A. C., Bouwman, J., Beckwith, M., van Dishoeck, E. F., Oberg, K. L, 
Pontoppidan, K. M., Linnartz, H., Blake, G. A., Evans, N. J., & Lahuis, F. 2010, ApJ, 718, 1100 

Bourke, T. L., Hyland, A. R., Robinson, G., James, S. D., & Wright, C. M. 1995, MNRAS, 276, 1067 



38 



Bourke, T. L., Myers, P. C, Evans, II, N. J., Dunham, M. M., Kauffmann, J., Shirley, Y. L., Crapsi, A., 
Young, C. H., Huard, T. L., Brooke, T. Y., Chapman, N., Cieza, L., Lee, C. W., Teuben, P., & 
Wahhaj, Z. 2006, ApJ, 649, L37 

Brinch, C, Crapsi, A., J0rgensen, J. K., Hogerheijde, M. R., & Hih, T. 2007, A&A, 475, 915 

Burkert, A., & Bodenheimer, P. 1993, MNRAS, 264, 798 

Casehi, P., Benson, P. J., Myers, P. C, & Tafaha, M. 2002, ApJ, 572, 238 

Cassen, P., & Moosman, A. 1981, Icarus, 48, 353 

Chapman, N. L., & Mundy, L. G. 2009, ApJ, 699, 1866 

Chen, X., Launhardt, R., & Henning, T. 2007, ApJ, 669, 1058 

Chiang, H., Looney, L. W., Tobin, J. J., & Hartmann, L. 2010, ApJ, 709, 470 

Chini, R., Ward-Thompson, D., Kirk, J. M., Nielbock, M., Reipurth, B., & Sievers, A. 2001, A&A, 369, 155 

Crapsi, A., Casehi, P., Walmsley, C. M., Myers, P. C, Tafaha, M., Lee, C. W., & Bourke, T. L. 2005, ApJ, 
619, 379 

Crapsi, A., Casehi, P., Walmsley, C. M., Tafaha, M., Lee, C. W., Bourke, T. L., & Myers, P. C. 2004, A&A, 
420, 957 

Danby, C, Flower, D. R., Valiron, P., Schilke, P., & Walmsley, C. M. 1988, MNRAS, 235, 229 

Di Francesco, J., Myers, P. C, Wilner, D. J., Ohashi, N., & Mardones, D. 2001, ApJ, 562, 770 

Djupvik, A. A., Andre, P., Bontemps, S., Motte, F., Olofsson, G., Galfalk, M., & Floren, H. 2006, A&A, 

458, 789 

Dunham, M. M., Evans, N. J., Terebey, S., Dullemond, C. P., & Young, C. H. 2010, ApJ, 710, 470 

Dzib, S., Loinard, L., Mioduszewski, A. J., Boden, A. F., Rodriguez, L. F., & Torres, R. M. 2010, ApJ, 718, 
610 

Enoch, M. L., Evans, N. J., Sargent, A. L, & Glenn, J. 2009, ApJ, 692, 973 

Enoch, M. L., Glenn, J., Evans, II, N. J., Sargent, A. L, Young, K. E., & Huard, T. L. 2007, ApJ, 666, 982 

Evans, N. J., Dunham, M. M., J0rgensen, J. K., Enoch, M. L., Merin, B., van Dishoeck, E. F., Alcala, J. M., 
Myers, P. C, Stapelfeldt, K. R., Huard, T. L., Allen, L. E., Harvey, P. M., van Kempen, T., Blake, 
G. A., Koerner, D. W., Mundy, L. G., Padgett, D. L., & Sargent, A. I. 2009, ApJS, 181, 321 

Froebrich, D. 2005, ApJS, 156, 169 

Fuller, G. A., Lada, E. A., Masson, C. R., & Myers, P. C. 1995, ApJ, 453, 754 

Fuller, G. A., & Wootten, A. 2000, ApJ, 534, 854 

Goldsmith, P. F., & Langer, W. D. 1999, ApJ, 517, 209 

Goodman, A. A., Benson, P. J., Fuller, G. A., & Myers, P. C. 1993, ApJ, 406, 528 

Gramajo, L. V., Whitney, B. A., Gomez, M., & Robitaille, T. P. 2010, AJ, 139, 2504 

Gueth, F., & Guilloteau, S. 1999, A&A, 343, 571 

Gueth, F., Guilloteau, S., & Bachiller, R. 1996, A&A, 307, 891 



39 



Haisch, Jr., K. E., Lada, E. A., & Lada, C. J. 2001, ApJ, 553, L153 

Hartmann, L., Calvet, N., Allen, L., Chen, H., & Jayawardhana, R. 1999, AJ, 118, 1784 

Harvey, P. M., Chapman, N., Lai, S., Evans, II, N. J., Allen, L. E., J0rgensen, J. K., Mundy, L. C, Huard, 
T. L., Porras, A., Cieza, L., Myers, P. C, Merin, B., van Dishoeck, E. F., Young, K. E., Spiesman, 
W., Blake, G. A., Koerner, D. W., Padgett, D. L., Sargent, A. L, & Stapelfeldt, K. R. 2006, ApJ, 
644, 307 

Heitsch, F., Slyz, A. D., Devriendt, J. E. C, Hartmann, L. W., & Burkert, A. 2006, ApJ, 648, 1052 

Herbig, G. H., & Jones, B. F. 1983, AJ, 88, 1040 

Hernandez, J., Hartmann, L., Megeath, T., Gutermuth, R., Muzerolle, J., Calvet, N., Vivas, A. K., Briceho, 
C, Allen, L., Stauffer, J., Young, E., & Fazio, G. 2007, ApJ, 662, 1067 

Hirota, T., Honma, M., Imai, H., Sunada, K., Ueno, Y., Kobayashi, H., & Kawaguchi, N. 2011, PAS J, 63, 1 

Ho, P. T. P., & Townes, C. H. 1983, ARA&A, 21, 239 

Hogerheijde, M. R., van Dishoeck, E. F., Blake, G. A., & van Langevelde, H. J. 1998, ApJ, 502, 315 

Johnstone, D., Rosolowsky, E., Tafaha, M., & Kirk, H. 2010, ApJ, 711, 655 

J0rgensen, J. K. 2004, A&A, 424, 589 

J0rgensen, J. K., Bourke, T. L., Myers, P. C, Di Francesco, J., van Dishoeck, E. F., Lee, C, Ohashi, N., 
Schoier, F. L., Takakuwa, S., Wilner, D. J., & Zhang, Q. 2007, ApJ, 659, 479 

J0rgensen, J. K., Harvey, P. M., Evans, II, N. J., Huard, T. L., Allen, L. E., Porras, A., Blake, G. A., Bourke, 
T. L., Chapman, N., Cieza, L., Koerner, D. W., Lai, S., Mundy, L. G., Myers, P. C, Padgett, D. L., 
Rebull, L., Sargent, A. L, Spiesman, W., Stapelfeldt, K. R., van Dishoeck, E. F., Wahhaj, Z., & 
Young, K. E. 2006, ApJ, 645, 1246 

Kauffmann, J., Bertoldi, F., Bourke, T. L., Evans, II, N. J., & Lee, C. W. 2008, A&A, 487, 993 

Keto, E., & Rybicki, G. 2010, ApJ, 716, 1315 

Kirk, J. M., Ward-Thompson, D., Di Francesco, J., Bourke, T. L., Evans, N. J., Merm, B., Allen, L. E., 
Cieza, L. A., Dunham, M. M., Harvey, P., Huard, T., J0rgensen, J. K., Miller, J. F., Noriega-Crespo, 
A., Peterson, D., Ray, T. P., & Rebuh, L. M. 2009, ApJS, 185, 198 

Koda, J., Sawada, T., Wright, M. C. H., Teuben, P., Corder, S. A., Patience, J., Scoville, N., Donovan 
Meyer, J., & Egusa, F. 2011, ApJS, 193, 19 

Lada, C. J. 1987, in lAU Symp. 115: Star Forming Regions, ed. M. Peimbert & J. Jugaku, 1-17 

Larson, R. B. 1969, MNRAS, 145, 271 

Launhardt, R., Nutter, D., Ward-Thompson, D., Bourke, T. L., Henning, T., Khanzadyan, T., Schmalzl, 
M., Wolf, S., & Zylka, R. 2010, ApJS, 188, 139 

Launhardt, R., Sargent, A., & Zinnecker, H. 2001, in Astronomical Society of the Pacific Conference Series, 
Vol. 235, Science with the Atacama Large Millimeter Array, ed. A. Wootten, 134 — h 

Lee, C, Hirano, N., Palau, A., Ho, P. T. P., Bourke, T. L., Zhang, Q., & Shang, H. 2009, ApJ, 699, 1584 

Lee, J., Bergin, E. A., & Evans, II, N. J. 2004, ApJ, 617, 360 



40 



Loinard, L., Torres, R. M., Mioduszewski, A. J., & Rodriguez, L. F. 2008, ApJ, 675, L29 

Loinard, L., Torres, R. M., Mioduszewski, A. J., Rodriguez, L. F., Gonzalez-Lopezlira, R. A., Lachaume, 
R., Vazquez, V., & Gonzalez, E. 2007, ApJ, 671, 546 

Looney, L. W., Tobin, J. J., & Kwon, W. 2007, ApJ, 670, L131 

Lovas, F. J. 1992, Journal of Physical and Chemical Reference Data, 21, 181 

Markwardt, C. B. 2009, in Astronomical Society of the Pacific Conference Series, Vol. 411, Astronomical 
Data Analysis Software and Systems XVIII, ed. D. A. Bohlender, D. Durand, & P. Dowler, 251 — h 

McCaughrean, M. J., Rayner, J. T., & Zinnecker, H. 1994, ApJ, 436, L189 

McKee, C. F., & Ostriker, E. C. 2007, ARA&A, 45, 565 

Menten, K. M., Reid, M. J., Forbrich, J., & Brunthaler, A. 2007, A&A, 474, 515 

Myers, P. C, Bachiller, R., Caselh, P., Fuller, G. A., Mardones, D., Tafalla, M., & Wilner, D. J. 1995, ApJ, 
449, L65+ 

Myers, P. C, Fuller, G. A., Goodman, A. A., & Benson, P. J. 1991, ApJ, 376, 561 

Pineda, J. E., Goodman, A. A., Arce, H. G., Caselh, P., Foster, J. B., Myers, P. C, & Rosolowsky, E. W. 
2010, ApJ, 712, L116 

Pudritz, R. E., & Norman, C. A. 1983, ApJ, 274, 677 

Raghavan, D., McAlister, H. A., Henry, T. J., Latham, D. W., Marcy, G. W., Mason, B. D., Gies, D. R., 
White, R. J., & ten Brummelaar, T. A. 2010, ApJS, 190, 1 

Rebuh, L. M., Stapelfeldt, K. R., Evans, II, N. J., J0rgensen, J. K., Harvey, P. M., Brooke, T. Y., Bourke, 
T. L., Padgett, D. L., Chapman, N. L., Lai, S., Spiesman, W. J., Noriega- Crespo, A., Merin, B., 
Huard, T., Allen, L. E., Blake, G. A., Jarrett, T., Koerner, D. W., Mundy, L. G., Myers, P. C, 
Sargent, A. L, van Dishoeck, E. F., Wahhaj, Z., & Young, K. E. 2007, ApJS, 171, 447 

Reipurth, B., Raga, A. C, & Heathcote, S. 1996, A&A, 311, 989 

Rich, J. W., de Blok, W. J. G., Cornweh, T. J., Brinks, E., Walter, F., Bagetakos, L, & Kennicutt, R. C. 
2008, AJ, 136, 2897 

Rodriguez, L. F., Reipurth, B., Raga, A. C, & Canto, J. 1998, Rev. Mexicana Astron. Astrofis., 34, 69 

Rydbeck, O. E. H., Sume, A., Hjalmarson, A., Ellder, J., Ronnang, B. O., & Kollberg, E. 1977, ApJ, 215, 
L35 

Sault, R. J., Teuben, P. J., & Wright, M. C. H. 1995, in Astronomical Society of the Pacific Conference 
Series, Vol. 77, Astronomical Data Analysis Software and Systems IV, ed. R. A. Shaw, H. E. Payne, 
& J. J. E. Hayes, 433-+ 

Schoier, F. L., van der Tak, F. F. S., van Dishoeck, E. F., & Black, J. H. 2005, A&A, 432, 369 

Scholz, A., Wood, K., Wilner, D., Jayawardhana, R., Delorme, P., Caratti O Garatti, A., Ivanov, V. D., 
Saviane, I., & Whitney, B. 2010, MNRAS, 409, 1557 

Scale, J. P., & Looney, L. W. 2008, ApJ, 675, 427 

Shirley, Y. L., Evans, II, N. J., Rawlings, J. M. C, & Gregersen, E. M. 2000, ApJS, 131, 249 



41 



Shu, F., Najita, J., Ostriker, E., Wilkin, F., Ruden, S., & Lizano, S. 1994, ApJ, 429, 781 

Shu, F. H. 1977, ApJ, 214, 488 

Shu, F. H., Adams, F. C, & Lizano, S. 1987, ARA&A, 25, 23 

Smith, R. J., Glover, S. C. O., Bonneh, I. A., Clark, P. C, & Klessen, R. S. 2011, MNRAS, 411, 1354 

Stanke, T., Smith, M. D., Gredel, R., & Khanzadyan, T. 2006, A&A, 447, 609 

Stutz, A. M., Rieke, G. H., Bieging, J. H., Balog, Z., Heitsch, F., Kang, M., Peters, W. L., Shirley, Y. L., & 
Werner, M. W. 2009, ApJ, 707, 137 

Tafalla, M., Mardones, D., Myers, P. C., CaselU, P., Bachiller, R., & Benson, P. J. 1998, ApJ, 504, 900 

Tafalla, M., Myers, P. C., Caselh, P., & Walmsley, C. M. 2004, A&A, 416, 191 

Tafalla, M., Myers, P. C., Mardones, D., & Bachiller, R. 2000, A&A, 359, 967 

Tanner, J. D., & Arce, H. G. 2011, ApJ, 726, 40 

Terebey, S., Fich, M., Noriega- Crespo, A., Padgett, D. L., Fukagawa, M., Audard, M., Brooke, T., Carey, 
S., Evans, N. J., Guedel, M., Hines, D., Huard, T., Knapp, G. R., McCabe, C, Menard, F., Monin, 
J., & Rebuh, L. 2009, ApJ, 696, 1918 

Terebey, S., Shu, F. H., & Cassen, P. 1984, ApJ, 286, 529 

Thompson, A. R., Moran, J. M., & Swenson, Jr., G. W. 2001, Interferometry and Synthesis in Radio 
Astronomy, 2nd Edition, ed. Thompson, A. R., Moran, J. M., & Swenson, G. W., Jr. 

Tobin, J. J., Hartmann, L., Calvet, N., & D'Alessio, P. 2008, ApJ, 679, 1364 

Tobin, J. J., Hartmann, L., & Loinard, L. 2010a, ApJ, 722, L12 

Tobin, J. J., Hartmann, L., Looney, L. W., & Chiang, H. 2010b, ApJ, 712, 1010 

Tsitah, A. E., Bourke, T. L., Peterson, D. E., Myers, P. C, Dunham, M. M., Evans, N. J., & Huard, T. L. 
2010, ApJ, 725, 2461 

Visser, A. E., Richer, J. S., & Chandler, C. J. 2002, AJ, 124, 2756 

Volgenau, N. H., Mundy, L. G., Looney, L. W., & Welch, W. J. 2006, ApJ, 651, 301 

Whitney, B. A., Wood, K., Bjorkman, J. E., & Cohen, M. 2003, ApJ, 598, 1079 

Young, K. E., Enoch, M. L., Evans, H, N. J., Glenn, J., Sargent, A., Huard, T. L., Aguirre, J., Golwala, S., 
Haig, D., Harvey, P., Laurent, G., Mauskopf, P., & Sayers, J. 2006, ApJ, 644, 326 

Yun, J. L. 1996, AJ, 111, 930 

Zhou, S., Evans, H, N. J., Koempe, C, & Walmsley, C. M. 1993, ApJ, 404, 232 



This preprint was prepared with the A AS L^T^]X macros v5.2. 



42 



NgH^ (1-0) IRAM 30m Intensity K-km/s 
+ IRAC 8 micron image 

fri ii ' ' u 



Centroid (km/s) 
2.55 2.6 2.65 2.7 2.75 



FWHM Line Width (km/s) 
0.2 0.3 0.4 0.5 0.6 0.7 0.£ 




-< °r 



-50 ^ -100 100 / 50 



2.55 2.6 2.65 2.7 2.75 



IRAC 8 micron image \ 




^H 




1 


- 


1 






ii 


In 


/ 


i. 


-Ti 


^ 


^ 




- i 


w 


0.05 pc 


- 




1 


, 1 , 


• - 

1 1 




Aa ["] 



Aa [■'] 



Aa ["] 



Fig. 1. — LI 157- The toy row shows the data from the IRAM 30m and the bottom row shows 
the PdBI data. In the left column^ the 8/im IRAC images with N2H+ (J = 1 ^ 0) integrated 
intensity contours overlaid. The IRAM 30m contours start at the lOa level and increase in 20a 
increments while the PdBI data start at the zbScr level and increase in 20a increments, see Tables 
6 - 8 for values of aj. The middle column shows the line-center velocity fit of the N2H+ emission 
across the envelope and the right column shows the FWHM linewidth. The red and blue lines 
mark the central axis of the outflow and their respective orientation in the plane of the sky; the 
X marks the position of the protostar from dust continuum and/or 24/im emission and the black 
arrows indicate the direction of the velocity gradient derived from the velocity fleld. The velocity 
gradient direction for the PdBI data could not be flt for L1157. The single-dish data show an 
ordered velocity gradient along the fllamentary envelope, following the high density region as it 
curves south. The interferometer map also reflects this velocity gradient, but there is a reversal 
with red-shifted emission just east of the protostar. There is also a red-shifted feature that appears 
in the interferometer and single-dish map just southeast of the protostar, this appears to be due to 
outflow interaction given its location along the cavity wall. Both the single-dish and interferometer 
maps have large linewidths near the protostar, which appears to be due to outflow interaction. 
Reference positions for the observations are listed in Tables 2 and 5. 
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Fig. 2. — L1157- The to'p and middle rows show the VLA NH3 (1,1) data while the bottom row 
shows the CARMA N2H+ (J = 1 ^ 0) data. In the left column, the 8/im IRAC images with 
integrated intensity contours overlaid. The VLA contours start at the zh3a level and increase in 
Qa increments while the CARMA data start at the ±3cr level and increase in 3a increments, see 
Tables 6 - 8 for values of aj. The middle column shows the line-center velocity of the gas across 
the envelope and the right column shows the FWHM linewidth. The red and blue lines mark the 
central axis of the outflow and their respective orientation in the plane of the sky; the X marks the 
position of the protostar from dust continuum and 24/im emission and the black arrows indicate the 
direction of the velocity gradient derived from the velocity field. Notice the close correspondence 
of the intensity distribution and line kinematics between the NH3 and N2H^ maps indicating that 
they are likely tracing the same material. Reference positions for the observations are listed in 
Tables 2 - 5. 
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Fig. 3. — L1157- Maps of N2H+ emission from the PdBI in L1157 in three velocity ranges; the 
blue and red lines denote the blue and red-shifted sides of the outflow ( J0rgensen et al.pOOT ). The 
line-center emission from 2.33-3.1 km s~^ is plotted as grayscale with light gray contours in units 
of 15, 30, 25, 60, 75a where a=0.027 K km s~^. The red contours are emission between 3.2 and 4 
km s~^ and the blue contours are emission between 1.57 and 2.2 km s~^and plotted in units of 3, 
6, 9, 12, 15, 21, 27a. The red- and blue-shifted emission near the protostar is clearly shifted along 
the blue-shifted side of the outflow and both the blue and red-shifted emission seem to outline 
the outflow cavity wall. The spatial shifts along the outflow in the blue and red-shifted emission 
strongly suggest that the outflow is entraining material from the inner envelope. 
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Fig. 4.— L1165- Same as Figure [T] but with CARMA data in the hoUom row. The IRAM 30m 
contours start at 5cr with 5a intervals; the CARMA data start at ±3cr with Sa intervals. The 
single-dish data show that the cloud has a fairly constant velocity away from the protostar, but 
near the protostar the kinematic structure is distinct; however, the single-dish linewidth map shows 
no indication of enhancement near the protostar. The velocity field from the interferometer map 
shows considerable detail, with a velocity gradient across the protostar nearly normal to the outflow. 
There is also enhanced linewidth near the protostar in the central envelope. Reference positions 
for the observations are listed in Tables 2-5. 
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Fig. 5. — L1165- The left panel shows the IRAC 8/xm image with CARMA HCO+ blue and red- 
shifted emission, summed over -3.5 to -2 km s~^ and -1.0 to km s~^, plotted as blue and red 
contours respectively. The contours levels are ±3, 6, and 8.25cr (a=0.175 K) for the blue-shifted 
emission and ±3, 9, 18, and 27a ((j=0.212 K) for the red-shifted emission. The blue and red-shifted 
emission from HCO^ is located symmetrically about the protostar, normal to the outflow. The 
dashed lines mark the regions where the position-velocity cut was taken and point to respective 
ends of the PV plot in the right panel. The position of the protostar/continuum source is marked 
with a white cross. The posit ion- velocity cut shows that the blue and red-shifted emission traces 
higher velocity material and there is a slight gradient of material going to higher velocity closer 
to the continuum source. The PV plot contours start at 3a and increase in 3a intervals (cr=0.2). 
The solid-blue curve represents Keplerian rotation (or infall) for a 0.5 Mq central object and the 
dashed-blue curve is for a 2.0 Mq central object. 
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NgH (1-0) IRAM 30m Intensity K-km/s 
+ IRAC 8 micron image 
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Fig. 6. — CB230- Same as Figure [11 but the interferometer data are VLA NH3 (1,1) observations. 
The IRAM 30m contours start at 3a with 10a intervals; the VLA data start at zbSa with 3a 
intervals. The line-center velocity from the single-dish data traces a fairly smooth velocity gradient 
across the envelope with enhanced linewidth near the protostar. The NH3 intensity is less extended 
directly east of the protostar associated with the cutoff of Sfim extinction, the rest of the envelope 
appears fairly flat. The line-center velocity from the NH3 emission traces a velocity gradient similar 
to the single-dish data; however, the shift from red to blue-shifted emission is quite abrupt and the 
transition region itself is curved. The NH3 line- width does not show much detail other than having 
its peak coincident with the highest intensity NH3 emission. 
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NgH"^ (1-0) IRAM 30m Intensity K-km/s 
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Fig. 7. — HH108IRS- Same as Figure [4} but zooming in on the single-dish data in the middle row 
and then the interferometer data in the bottom row. The IRAM 30m contours start at 3a with lOa 
intervals; the CARMA data start at ±3cr with 6a intervals. The N2H+ peak in this object appears 
offset from the protostar in the single-dish and interferometer map. The interferometer map further 
appears to be elongated normal to the outflow. The single-dish line-center velocity map appears to 
show a slight gradient normal to the outflow, but the interferometric velocity fleld starts red-shifted, 
becomes blue-shifted and then goes back to red-shifted. The single-dish linewidth maps show a 
slight enhancement near the protostar; the interferometer data on the other hand show a highly 
increased linewidth near the protostar as compared to the surrounding region. 
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NgH"^ (1-0) IRAM 30m Intensity K-km/s 
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Fig. 8.— HH108MMS- Same as Figure [7| The IRAM 30m contours start at 3a with lOa intervals; 
the CARMA data start at zb3cr with 3a intervals. The single-dish integrated intensity map shows 
HH108MMS as an extension along a large-scale filamentary structure that connects with HH108IRS. 
However, the interferometer map finds the N2H+ peak directly coincident with the protostar. The 
single-dish velocity field is fairly constant throughout the region of the envelope, but the large scale 
gradient is normal to its outflow. The CARMA velocity map then shows that there is a gradient 
along the outflow and the linewidth map shows increased linewidth along the outflow, this detail 
was absent in the single-dish map. 
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NgH"^ (1-0) IRAM 30m Intensity K-km/s 
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Fig. 9. — Serpens MMS3- Same as Figure[T]but with VLA NH3 (1,1) data and an additional zoom- 
in of the NH3 data in the bottom row. The IRAM 30m contours start at 5cr with lOa intervals; 
the VLA data start at ±3a with 5cr intervals. The N2H+ peak in the single-dish data is not on 
Serpens MMS3 but rather near a small cluster of young stars to the southwest, but there is an 
extension toward MMS3. The single-dish velocity field near MMS3 shows a red-shifted pocket 
associated with its position with a gradient extending along the direction of the large-scale filament 
shown in the top panels. The VLA NH3 map reveals more detail as the emission closely follows 
the filamentary structures in the region. The NH3 velocity map also shows the deep red-shifted 
pocket of emission next to the protostar and reveals more detail in the gradient along the filament. 
There is a corresponding increase in linewidth to the east of the protostar at the location of the 
red-shifted emission. 
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Fig. 10. — HH211- Same as Figure [4j The IRAM 30m contours start at 3cr with lOcr intervals; 
the CARMA data start at ±3cr with 3cr intervals. As shown in the left panels, the N2H+ emission 
mapped by the 30m and CARMA correlates very well with the 8/im extinction in the region. We 
also notice that the N2H+ peak is offset from the protostar to the southwest in the single-dish and 
interferometer data. The velocity field shows a linear gradient normal to the outflow of HH211; 
however, south of the protostar there is another N2H+ velocity component blue-shifted from the 



rest of the gas in the region. [Tanner fc Arce (2011) referred to this as the southwest extension. 
This transition region appears as artificially large line-width in the top right panel. The velocity 
field in the CARMA data also shows the linear gradient with a slight increase in linewidth near the 
protostar. 
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Fig. 11.— IRAS 16253-2429- Same as Figureji} Both the IRAM 30m and CARMA contours start at 
3(7 with 3a intervals. The single-dish N2H^ emission traces the roughly symmetric envelope of IRAS 
16253-2429 quite well with the peak emission just east of the protostar. The CARMA N2H+ map 
also traces the envelope well and the northern outflow cavity is prominent as an evacuated region, 
probably due to resolved-out emission. The single-dish velocity field shows a very small velocity 
gradient roughly normal to the outflow. The CARMA velocity map shows complex structure, 
but there is an overall gradient in the same direction as the single-dish map. In addition, there 
is a red-shifted feature at the same location of as the outflow cavity wall. The single-dish and 
interferometer linewidth maps show slight enhancements at the same locations, where the peak 
emission is present, but there also is a linewidth increase toward the edge of the envelope along the 
outflow in the single-dish map. Note that the linewidth in the envelope is extremely narrow, only 
^0.2 km s~^ over the regions where the outflow could not be interacting. 
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NgH"^ (1-0) IRAM 30m Intensity K-km/s 
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Fig. 12. — LI 152- Same as Figure [4j The IRAM 30m contours start at lOa with lOa intervals; the 
CARMA data start at ±3a with 3a intervals. The single-dish N2H+ emission traces two connected 
peaks over ^0.1 pc; one is star-less the other is adjacent to the protostar. The CARMA N2H+ 
map focuses on the protostellar clump and clearly shows that the protostar is offset from most of 
the N2H+ emission. Both velocity maps show a strong gradient normal to the outflow, while the 
star-less clump does not have much velocity structure. The linewidth maps show increases along 
the outflow; there is a feature in the single-dish map which is not exactly along the outflow but 
quite near it. The CARMA linewidth map shows increased linewidth that correlates very strongly 
with the jet-like emission shown in the 8/im image. 
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NgH"^ (1-0) IRAM 30m Intensity K-km/s 
+ IRAC 8 micron image 
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Fig. 13. — L1527- Same as Figure [4} but VLA NH3 data are also shown in the bottom row. The 
IRAM 30m contours start at 5a with 10a intervals; both the CARMA and VLA contours start at 
±3a with 3a intervals. We show the NH3 as well since they have slightly better resolution and 
because the N2H+ data only have moderate signal-to-noise. The single-dish N2H+ data are peaked 
north of the protostar, which is also the case for the CARMA data, but the VLA data are peaked 
south of the protostar. Both the VLA and CARMA data seem to trace a structure curving to the 
northeast which may be associated with the outflow cavity. The single-dish velocity fleld appears 
to have a component along the outflow and normal to it. The interferometer velocity maps reflect 
the single-dish velocities on the largest scales, but near the protostar there is a small-scale velocity 
gradient in the opposite direction of the large-scale gradient. There is little increase in linewidth 
in the single-dish data, but there appears to be increased linewidth near the protostar in both sets 
of interferometer data. 
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NgH"^ (1-0) IRAM 30m Intensity K-km/s 
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Centroid (km/s) 
10 10.1 10.2 10.3 10.4 10.5 10.6 



FWHM Line Width (km/s) 
0.3 0.4 0.5 0.6 0.7 O.E 





Aa [■'] \ / Aa ["] 

^ \ ^ Centroid (km/s) \ 

/ \ / 10 10.1 10.2 10.3 10.4 10.5 10.^ 

/NsH"^ (1-0) CARMA Intensity K-km/^ 
IRAC 8 micron image \ 



/ FWHM Line Width (km/s) \ 

0.2 0.3 0.4 0.5 0.6 , 





0.05 pc 



I I i 


Aa ["] 



/ 




0.05 pc 



Aa ["] 



Fig. 14. — RN043- Same as Figure [4j The IRAM 30m contours start at 3a with 10a intervals; the 
CARMA contours start at =b3cr with 3a intervals. The single-dish intensity map appears mostly 
unresolved, while the CARMA data trace the 8/im extinction structures closely with a depression 
of emission coincident with the protostar. The single-dish velocity field shows a gradient mostly 
normal to the outflow and the interferometer data reveal a sharp velocity shift between the east 
and west sides of the envelope. This appears as a line of broad linewidth (an artifact from fitting) 
where the two components are viewed on top of each other; this is associated with a region of large 
linewidth in the single-dish map. 
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Fig. 15.— RN043- Same as Figure[5| The contours levels are ±5, 8, 11, 14, and 17a {a = 0.12 K) 
for the red and blue-shifted emission. The emission morphology of the HCO^ in both the integrated 
intensity map on the left and in the PV plot are quite similar to LI 165. The total velocity extent 
of ^ ±2 km s~^ is similar in magnitude to LI 165. However, the red and blue peaks are much closer 
together due to the increased distance of RN043. The PV plot contours start at 3a and increase 
in 3a intervals (cr=0.2671). The solid-blue curve represents Keplerian rotation (or infall) for a 0.67 
Mq central object and the dashed-blue curve is for a 2.67 Mq central object. 
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Fig. 16. — IRAS 04325+2402- Same as Figure [TJ but we only have single-dish data for this object. 
The IRAM 30m contours start at lOa with lOcr intervals. The N2H+ contours correlate very 
well with the Sfim extinction and the peak is actually on an apparent star-less core north of the 
protostar. There is a tail of emission wrapping towards the protostar, but at lower intensity levels. 
The velocity field shows a gradient roughly normal to the outflow and there appears to be smaller- 
scale structure in the velocities near the protostar. The linewidth map does not show a large 
increase near the protostar but shows two peaks that are located along the outflow indicating a 
possible interaction. 
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NgH"^ (1-0) IRAM 30m Intensity K-km/s 
+ IRAC 8 micron image 



Centroid (km/s) 
5.2 5.3 5.4 5.5 5.6 5.7 



FWHM Line Width (km/s) 
0.2 0.3 0.4 0.5 0.6 0.7 0.8 



- * f M 


^^, , , , 1 , , ,_ 


\*\ ( 1 / // 


sS^^t^ 




^f'F 




oT?^p\» *^^^ \ 


'^^/l 1 1 \ \ IM 1 1 Wl 1 I \ 1 \l 



m> > 


V 


^^^^ 1 1 ^ ^^Y y 




' -^^7 


-<^t 


^B/ pc / 


1 " 
) 1 J 
^^^-^ - 


^^^/i 1 1 1 1 1 1 1 'i 1 w r 1 1 \ 1 



Aa [■'] 



50 \ / 50 

\ / 

/ \ / 

/ NHg (1,1) VLA Intensity K-km/s ^ 
+ IRAC 8 micron image ^ 







Centroid (km/s) 
5.2 5.3 5.4 5.5 5.6 5.7 




FWHM Line Width (km/s) \ 

0.2 0.3 0.4 0.5 0.6 0.7 0.8 * 



.^ 




Fig. 17. — L483- Same as Figured) using archival NH3 data from Fuller h Wootten (2000). The 



IRAM 30m contours start at 20cr with 20cr intervals; the VLA contours start at ±3cr with Go- 
intervals. The N2H^ and NH3 emission trace the envelope seen in 8/im extinction quite closely. 
The velocity field shows a gradient that is at an angle 45° from the outflow in the single-dish map. 
The VLA velocity map shows a similar large-scale feature but there is a pocket of blue-shifted 
gas coincident with the protostar and there is a red-shifted pocket directly north of the protostar. 
The single-dish linewidth map shows increased linewidth at the location of the most rapid velocity 
changes; the VLA linewidth map only shows a couple of features where there is a sharp transition 
in the velocity components. This is due to only having 0.3 km s~^ channel width in the NH3 data. 
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NgH (1-0) IRAM 30m Intensity K-km/s 
+ IRAC 8 micron image 
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Fig. 18. — L673- Same as Figure [6] with an additional zoom-in on the NH3 data in the bottom row. 
The IRAM 30m contours start at 20a with 20cr intervals; the VLA contours start at zbScr with Go- 



intervals. Though only three protostars are marked; there are likely seven in the region (Tsitali 



et al.|[20Tol ); the southern- most protostar marked is actually composed of three sources in Ks-band 
imaging (Tobin et al. in prep.). The single-dish N2H+ emission fills the map, but there is increased 
emission where the S/xm extinction is most prominent. The NH3 map only picks up the densest 
regions around the protostars. There is a gradient in the velocity map along the filament, though 
there appears to be features associated with the two marked protostars in both the single-dish and 
VLA maps. Both velocity maps also show a second velocity component to the north. The linewidth 
maps show increased linewidth at the transition region between the two components and the VLA 
map shows increased linewidth near the protostars. 
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NgH"^ (1-0) IRAM 30m Intensity K-km/s 
+ IRAC 8 micron image 
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Fig. 19. — L1521F- Same as Figure [6j The IRAM 30m contours start at 3cr with 3cr intervals; 
the VLA contours start at dz3cr with 3cr intervals. The single-dish data trace the region of 8/xm 
extinction well and the NH3 data trace the inner, high density regions. The velocity structure 
from the single-dish map shows weak evidence for a gradient normal to the outflow. The NH3 
data do show a gradient, but the emission near the protostar is all red-shifted and the blue-shifted 
emission is far (10^^ 1400 AU) from the protostar. However, as noted in the text, the NH3 emission 
is optically thick and may not fully trace the inner envelope kinematics. The single-dish linewidth 
map has very little structure while the NH3 linewidths are peaked coincident with the blue-shifted 
area of the velocity field. 
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NgH"" (J=l-0) CARMA Intensity K-km/s 
+ IRAC 8 micron image 



Centroid (km/s) 
7.9 8 8.1 8.2 8.3 




NHgD (1(1,1)0-1(0,1)0) CARMA Intensity K-km/s 
+ IRAC 8 micron image 




FWHM Line Width (km/s) 
0.2 0.3 0.4 0.5 0.6 0.7 0.8 






li 


- 


II II 

0.05 pc 




:'^ 


^ 




• 



Centroid (km/s) 
7.9 8 8.1 8.2 8.3 


^■^ 




P 




0.05 pc 




, 1 , 


, 1 , 


1 1 1 1 1 1 1 1 1 

•- 



20 -20 

Aa ["] 

FWHM Line Width (km/s) 
0.2 0.3 0.4 0.5 0.6 0.7 0.£ 







II 








: 3 


0.05 pc 


-^ 


^ 


1 


1 1 1 1 


1 1 1 

• - 



La ["] 



Aa ["] 



Aa ["] 



Fig. 20. — Perseus 5- The toy row shows the CARMA N2H+ ( J = 1 ^ 0) data and the bottom row 
shows the NH2D (li^i -^ lo,i) data also from CARMA; the contours start at zbScr with Qa intervals 
for N2H+ and ±3cr with 3a intervals for NH2D. Both molecules have approximately the same noise 
level. The integrated intensity data for both molecules are overlaid on the IRAC Sfim images. 
NH2D appears to trace the S/xm extinction best while N2H+ is centrally peaked. The line-center 
velocity of each molecular line is shown in the center column, N2H^ and NH2D indicate that there 
may be a gradient normal to the outflow, but that there is also a gradient in the direction of the 
outflow. The linewidth plots in the right panels are complex, all show an enhancement away from 
the protostar along the outflow; NH2D in particular shows an increase through the envelope near 
the outflow axis. 
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NgH"^ (1-0) IRAM 30m Intensity K-km/s 
+ IRAC 8 micron image 



Centroid (km/s) 
9 7 7.1 7.2 7.3 7.4 



FWHM Line Width (km/s) 
0.2 0.3 0.4 0.5 0.6 0.7 




/ 50 -50 \ / 50 -50 \ / 50 -50 ' 

I Aa ["] \ I Aa ["] \ / Aa ["] \ 

/ \ / Centroid (km/s) ^ / FWHM Line Width (km/s) 

\ 6.8 6.9 7 7.1 7.2 7.3 7.4 \ / 0.2 0.3 0.4 0.5 0.6 0.7 ^ 

' NH3 (1,1) VLA Intensity K-km/s y ' 
/ + IRAC 8 micron image ^ 



m 




^1 1 1 1 11 




1 1 1 1 


1 1 1 


' 1 


_ 


- 


^ 


''/ 


T 


- 


i^ 


Mt 






- 


1 1 


0.05^ 




• 

1 1 


- 


,1,1 


1 1 1 



Aa ["] 



Aa [■'] 



Aa ["] 



Fig. 21.— IRAS 03282+3035- Same as Figure JG] The IRAM 30m contours start at 3a with 10a 
intervals; the VLA contours start at ±3a with 5a intervals. The single-dish N2II+ emission closely 
traces the entire area viewed in S/xm extinction, peaking near the protostar. The VLA NH3 data 
are picking up the densest region near the protostar. The velocity gradient is mostly along the 
outflow in the single-dish and VLA NH3 velocity map. The linewidth in the single-dish data peaks 
to the west of the protostar, somewhat along the outflow. The NH3 linewidth map has peaks about 
the protostar, normal to the outflow, and directly north, while the rest of the map has a small 
linewidth. 
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Fig. 22. — Channel maps of NH3 (1,1) emission (contours) overlaid on CO ( J = 1 ^ 0) emission 



from Arce &; Sargent (2006), showing the relationship between NH3 emission and the outflow cavity. 
The top and bottom rows show the red- and blue-shifted component that are unblended; the NH3 
emission in the middle panels is blended. Contours start at ±3cr, and then increase in 3a intervals. 
The negative contours are plotted as dashed lines, reflecting the loss of large-scale structure. The 
unblended velocity channels show very narrow structures and in the bottom row, the blue-shifted 
components are almost north-south in orientation. Much of the NH3 emission south of the protostar 
appears unlikely to be influenced by the outflow. The position of the protostar is marked with the 
cross and the outflow axis is marked bye the blue and red lines denoting the outflow orientation 
within the plane of the sky. 
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NgH"" (1-0) IRAM 30m Intensity K-km/i 
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Fig. 23.— HH270 VLAl- Same as Figure [i) The IRAM 30m contours start at 3a with 6a intervals; 
the CARMA contours start at =t3cr with 3a intervals. The single-dish map shows strong emission 
coincident with the 8/xm extinction extended along the outflow. The CARMA N2H+ map shows a 
similar morphology, with no emission coincident with the protostar. The single-dish velocity fleld 
shows a gradient normal to and along the outflow; however, the CARMA velocity map shows that 
the gas velocity is relatively constant normal to the outflow, but along the outflow axis there is a 
clear gradient with more red-shifted emissions toward the outflow axis. The velocity contours trace 
a similar shape as the outflow cavity; the faint scattered light of the outflow cavity can be seen in 
the 8/im image. The single-dish linewidth map shows a peak along the outflow southwest of the 
protostar, this appears to be the unresolved velocity gradient along the outflow as the CARMA 
map does not have increased linewidth along at the same location. 
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Fig. 24. — Line center velocity cuts taken across each protostellar envelope, normal to the outflow. 
The diamonds are the single-dish data and the plus-signs are the interferometer data. Each point 
for the single-dish data is the average of the 27^' beam surrounding each point. Each interferometer 
point is the average of points within the semi-major axis of the synthesized beam. The velocity 
errors at each point are generally <0.05 km s~^, approximately the size of the symbols. The solid- 
lines are the linear fits to the single-dish velocity data between ±30^' for most sources. Serpens 
MMS3 was fit between 60'' and 0'' and L1152 was fit between 40'' and -20". 
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Fig;. 24 — cont'd. 
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Fig. 24 — cont'd. 
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Fig. 25. — Histogram plots of measured velocity gradients from single-dish data. The left panel 
shows the velocity gradients derived from one-dimensional cuts across the velocity field, taken nor- 
mal to the outflow. The middle panel shows the velocity gradients derived from a two-dimensional 
fit to the velocity field. The right panel shows the velocity gradients derived from the velocity dif- 
ference at ±10000 AU from the protostar, also normal to the outflow. The distributions from the 
one- and two-dimensional fits are comparable while the two point method is skewed toward smaller 
gradients; this difference is likely due to the gradients fits picking up on higher velocity emission 
that sometimes turns over toward lower velocities by ^ 10000 AU. The middle panel include the 
two-dimensional fit for L673 which is absent from the left and right panels. 
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Fig. 26. — Histogram plots of measured velocity gradients from interferometer data. The left 
panel shows the velocity gradients derived from one-dimensional cuts across the velocity field, taken 
normal to the outflow. The right panel shows the velocity gradients derived from a two-dimensional 
fit to the velocity field. The differences in these distributions likely result from the two-dimensional 
method having to fit all the data where complexities in the velocity field may reduce the gradient 
fit. The left panels include L1157 and Serpens MMS3 while they are absent from the right panels. 
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Fig. 27. — Histogram plots of velocity gradient position angle offset relative to outflow position angle 
for single-dish (left panel) and interferometric (right panel) velocity gradient measurements. An 
offset of 90° indicates that the gradient is normal to the outflow and most of the velocity gradients 
are within 45° of normal to the outflow (dashed line); however, there is substantial dispersion in 
this relationship. The shift in number toward 90° in the interferometer observations may reflect 
that smaller-scale motion is becoming more ordered. The left panels include LI 157 and Serpens 
MMS3 while they are absent from the right panels. 
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Fig. 28. — Histogram plot of single-dish velocity gradient position angle minus the interferometric 
velocity gradient position angle. Most velocity gradients at large and small-scales are within 45° of 
each other, indicating that the line of sight velocities at large-scales reflect similar velocity structure 
at small-scales. This plot does not include LI 157, Serpens MMS3, or L673. 
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Fig. 29. — Plot of escape velocity at 10000 AU versus velocity shifts measured at 10000 AU radii 
(diamonds) and the average N2H+ FWHM from the single-dish data (squares). The escape velocity 
is calculated using the total mass of the envelope from the 8/xm extinction data in Paper I plus 
a 1 Mq solar mass central object; the error bars in the calculated escape velocity reflect a 50% 
uncertainty in total mass. This shows that the envelopes are consistent with being gravitationally 
bound on large-scales and that the envelopes are not supported by rotation, turbulence, and/or 
thermal pressure. The protostar closest to rotational support is HH211; its substantial rotation 



was also highlighted by Tanner & Arce (2011). 
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Fig. 30. — Plot of single-dish N2H+ line profiles at three positions in HH211 showing the transi- 
tion between different velocity components. The middle spectrum is in the region where the two 
components are blended, demonstrating that the lines are not extremely wide. 
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Fig. 31. — IRAC 8/xm images of sources observed, with CARMA A=3mm continuum contours 
overlaid. For most protostars, the continuum emission is directly coincident with the 8/im point 
sources; 70/xm in the case of HH108IRS/MMS. The contours for Perseus 5, HH270, IRAS 16253- 
2429, HH108MMS, and L1152 are ±2.5cr, 3,4,5,. ..,9cr and then increase in 3a intervals. The contours 
for HH211, L1527, RN043, HH108IRS, and L1165 start at ±3a and increase by 3a. The contours 
for L1157 are ±3, 6, 12, 24, 36, 48, 60, 72, 84, and 96cr. 
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